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ABSTRACT 


This  thesis  evaluates  the  ability  of  the  IRTOOL  computer  simulation  program  to 
predict  mirages.  Using  identical  input  conditions  taken  from  the  MAPTIP  experiment 
database,  predicted  Minimum  Mirage  Range  (MMR)  and  Maximum  Intervision 
Range  (MFVR)  from  both  the  IRTOOL  and  IRBLEM  models  were  extracted  and 
compared  with  the  measurements  recorded  in  the  database.  By  comparison  of  the 
algorithms  it  was  found  that  discrepancies  in  IRTOOL  mirage  prediction  could  be 
ascribed  to  the  input  function  for  significant  ocean  wave  height,  which  gave  values 
much  greater  than  measured  or  used  in  IRBLEM.  For  a  significant  wave  height  close 
to  the  measured  value  the  IRTOOL  predictions  were  in  very  close  agreement  with 
observation  and  with  IRBLEM.  IRTOOL  predictions  were  in  all  cases  within  2.7  km 
and  in  most  cases  within  1.3  km  of  the  measurements  for  all  ranges  varying  from 
about  7-26  km.  The  strong  temperature  gradient  predicted  by  the  model  within  a  few 
meters  of  the  water  surface,  uncertainties  in  the  measured  range,  and  the  variation  of 
0.8  to  2°C  in  Air  Sea  Temperature  Difference  are  sufficient  to  account  for  the 
observed  deviations.  Differences  between  the  model  predictions  and  some  of  the 
problems  encountered  are  also  discussed. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Modem  sensor  and  weapon  systems  rely  on  propagation  of  electromagnetic  or 
electro-optic  energy  in  an  inhomogeneous  atmosphere.  In  the  coastal  environment, 
complex  spatial  structures  cause  a  great  variability  in  the  atmospheric  parameters  that 
affect  the  detection  of  low  flying  targets.  One  of  the  parameters  among  others  which  is 
responsible  for  the  poor  performance  of  the  electro-optical  systems  is  the  gradient  of 
refractivity  caused  by  large  vertical  changes  in  temperature  and  humidity,  especially  near 
the  ocean  surface.  Shipboard  optical  systems  are  used  currently  for  the  detection  and 
acquisition  of  low  flying  air  targets.  The  range  performance  of  IRST  (Infrared  Search  and 
Track)  sensors  is  partly  determined  by  the  propagation  of  the  atmosphere.  This  is  why  the 
Navy  has  a  particular  interest  in  determining  the  electro-optical  properties  of  the  lower 
layers  of  the  marine  atmosphere.  In  recent  years  many  scientific  teams  from  different 
countries  have  been  trying  to  create  computer  based  simulation,  modeling,  performance 
prediction  systems  and  techniques,  taking  into  account  the  relevant  atmospheric 
parameters.  These  parameters  can  be  obtained  through  sensing,  numerical  modeling  or  a 
combination  of  both,  so  as  to  assess  refractive  effects  on  detection  range  of  low  level 
targets  including  mirage  and  multiple  image  generation.  In  spite  of  the  efforts  that  have 
been  put  forth  to  create  computer  models  which  give  results  very  close  to  the 
experimental  data,  the  detection  of  low  flying  anti-ship  missiles  (sea-skimmers)  with  high 
detection  probability  and  low  false  alarm  rates  is  still  a  challenge. 
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B.  STUDY  OBJECTIVES 

The  IRTOOL  (Infrared  tool)  computer  program  is  a  simulation  computation 
model  environment  joining  a  number  of  component  models  from  various  sources  into  a 
unified  package  that  produces  both  a  visual  simulation  and  engineering  calculations  for 
scenarios  appropriate  to  the  Infrared  Search  and  Track  systems  (IRST).  It  has  been 
developed  by  Arete  Corporation  and  NSWC  (Naval  Surface  Warfare  Center)  under  the 
sponsorship  of  the  Office  of  Naval  Research  in  support  of  the  Infrared  Analysis  Modeling 
and  Measurements  Program  (IRAMMP)  to  assist  in  the  design  phase  of  the  next  IRST 
generation.  It  is  a  DoD  generated  program  not  for  sale  and  commercial  use.  A  group  of 
scientists  from  the  Canadian  Defense  Research  Establishment  Valcartier  (DREV)  applied 
IRTOOL  to  data  they  obtained  off  the  Netherlands  coast  in  the  MAPTIP  measurements  of 
1993,  and  concluded  that  IRTOOL  did  not  predict  mirages  in  cases  in  which  they  were 
found  experimentally.  They  then  decided  to  forego  the  use  of  IRTOOL  completely  to 
concentrate  on  development  of  the  local  program  IRBLEM  (L.Forand  private 
communication) . 

The  IRBLEM  is  another  computer  code  (a  software  package)  for  the  computation 
of  atmospheric  effects  on  EO  band  systems  in  the  marine  surface  layer.  The  ray  refraction 
module  (L(W)WKD)  of  IRBLEM  has  been  evaluated  in  the  past  with  the  French 
CELAR’s  (Centre  d’Electronique  de  l’Armement)  PIRAM  (Profils  d’lndice  de  Refraction 
en  Atmosphere  Marine)  ray  refraction  model  using  the  database  of  the  MAPTIP 
experiment.  The  results  for  both  models  were  very  satisfactory. 


The  main  objective  of  this  research  is  the  evaluation  of  IRTOOL  for  visible/IR 
refraction  and  mirage  formation  in  the  marine  boundary  layer  as  well  as  the  comparison 
of  the  above  mentioned  code  with  IRBLEM.  For  the  comparison  between  those  two 
models  the  experimental  database  from  the  MAPTIP  trial,  which  took  place  in  the 
Netherlands  in  1993,  by  French,  Canadian  and  German,  teams,  was  used.  Using  the  same 
set  of  input  conditions  for  both  models  we  compare  the  predicted  results  with  the 
experimental  data. 

A  study  has  been  made  in  the  past  year  in  the  evaluation  and  comparison  of  the 
IRTOOL  refraction  model  with  the  above  mentioned  programs  IRBLEM  and  PIRAM 
[Ref.  2].  The  conclusion  of  this  effort  was  not  very  satisfactory  (i.e.  IRTOOL  did  not 
predict  mirages  in  cases  in  which  they  were  found  experimentally).  So,  the  second 
objective  of  this  thesis  is  the  use  of  as  many  cases  (data  sets)  as  possible,  so  as  to  cover 
all  the  possible  situations  and  seek  the  modification  of  the  IRTOOL  code,  if  that  is 
necessary,  which  may  lead  to  greatly  improved  performance  under  mirage  conditions. 
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II.  THEORY 


A.  BACKGROUND 

The  index  of  refraction  n,  is  a  dimensionless  parameter  given  by 

n=-  (1) 

v 

where  c  is  the  phase  velocity  of  the  electromagnetic  wave  (EM)  in  vacuum  given  by 

c=-rL==3xl08m/s  (2) 

4£oVo 

v  is  the  phase  velocity  of  the  EM  wave  in  a  homogeneous  medium  given  by 

(3) 

s0  =8.85  x  10  ~12  F/m  is  called  the  permittivity  of  free  space 

/j.0  =  1.26x  10"6  H/m  is  called  the  permeability  of  free  space  and 

s  and  n  are  the  permittivity  and  permeability  of  the  medium  respectively. 

Due  to  the  fact  that  in  general  /j.  =  /u0  and  s  >£0,  the  light  propagates  slowly  through 

matter.  The  EM  properties  of  the  medium  are  related  to  the  index  of  refraction  according 
to  Eqs.  1,  2,  3  by  the  equation: 

n= 


1 
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Any  change  to  propagation  velocity  due  to  the  interaction  with  the  medium  results  in  a 
bending,  called  refraction,  of  the  signal  path  (i.e.,  the  bending  of  the  rays).  The  index  of 
refraction  (n)  is  a  function  of  frequency  (color).  This  phenomenon  is  called  dispersion  and 
the  supporting  medium  is  called  dispersive.  The  dependence  of  n(co)  on  the  frequency  for 
different  optical  materials  for  the  visible  range  is  shown  in  Figure  1 . 


Figure  1 .  Dispersion  curves  for  several  different  optical  materials  commonly  used  for 
lenses  and  mirrors  in  the  visible  range  [Ref.l]. 


An  important  law,  named  after  the  Dutch  mathematician  Willebrord  van  Roijen 
Snell,  states  that  the  product  of  the  refractive  index  and  the  sine  of  the  angle  of  incidence 
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of  a  ray  in  one  medium  is  equal  to  the  product  of  the  refractive  index  and  the  sine  of  the 


angle  of  refraction  in  a  successive  medium  i.e.,  n,  sin#,  =n2  sin#2  (Figure  2).  So,  when  a 


ray  enters  a  medium  of  a  higher  refractive  index  it  will  bend  toward  the  normal  and  when 


a  ray  enters  a  medium  of  lower  refractive  index  it  will  deviate  away  from  the  normal 


(Figure  2).  The  denser  a  material  is,  the  higher  is  its  refractive  index.  Given  that  the  phase 


c 

velocity  of  the  EM  wave  is  v  =— ,  in  a  denser  medium  the  velocity  of  light  will  be  lower. 

n 


Figure  2.  Refracted  ray  passing  through  different  layered  media  of  increasing  refractive 
index  [Ref.  2], 
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As  a  result  we  expect  that  light  traveling  through  air  travels  more  slowly  as  the  density  of 
the  air  increases.  In  Figure  3  the  curve  is  the  ray  path  of  a  traveling  light  wave  and  the 
lines  perpendicular  to  this  curve  represent  wave  fronts.  The  curved  path  that  the  ray 
follows  is  due  to  the  decrease  of  the  air’s  density  with  height.  The  lower  part  of  the  wave 
front  is  lower  in  the  atmosphere  where  the  density  is  higher,  so  it  moves  more  slowly  than 
the  upper  part,  causing  the  ray’s  bending.  For  an  observer  on  earth  the  starlight  appears  to 
be  coming  from  a  higher  elevation  than  the  true  one  [Ref.  3].  This  is  analogous  to  the 
spoke  of  a  spinning  wheel  where  the  outer  part  of  the  wheel  is  moving  faster  than  the 
inner  part. 


In  order  to  calculate  the  ray  bending  we  will  use  the  model  of  Figure  4  [Ref.  4], 


Figure  4.  Beam  curvature  in  an  slightly  inhomogeneous  atmosphere  [Ref.  4]. 


The  adjacent  ray  paths  separated  by  dz  have  length  s  and  s+ds  and  also  the  refractive 
indices  of  the  air  along  the  paths  are  n  and  n+dn  (dn<0).  The  optical  path  that  each  ray 
travels  is  stationary  with  respect  to  variations  of  that  path  (Fermat’s  principle).  Hence, 

d(ns)  _  ds  dn 


•=n— +s— =0 
dz  dz  dz 


(5) 


ds 

But,  —  is  the  change  in  the  direction  of  the  wavefront  so  that  the  direction  change  per 

dz 


unit  path  length  (i.e.  the  path  curvature)  is: 


1  ds  l  dn 

cr= - = - [Ref.  4] 

s  dz  n  dz 


(6) 
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Because  —  is  negative  in  our  case,  the  ray  bends  downward.  This  is  usually  the  case  in 
dz 

the  atmosphere,  so  for  convenience  our  convention  will  be  that  such  a  downward  oriented 
curve  has  a  positive  radius.  The  value  of  the  refractive  index  of  air  is  very  close  to  unity, 
and  so  the  radius  of  curvature  of  the  beam  curvature  is  given  by: 


_  J__  1 

o"  ,dn 


(7) 


<*> 


Here,  it  is  very  important  to  note  that  the  ray  bending  has  a  radius  of  curvature  that 
depends  on  the  gradient  of  the  index  of  refraction  and  not  on  the  absolute  magnitude  of 
the  index  of  refraction  itself. 


B.  ELECTROMAGNETIC  SPECTRUM 

Many  times  in  this  thesis  we  will  refer  to  the  EM  spectrum  and  its  wavelength 
bands.  So,  it  will  be  very  useful  to  introduce  the  EM  spectrum.  The  EM  spectrum  is 
divided  into  six  generic  spectral  regions.  Four  of  these  regions  are  associated  with 
infrared  imaging  systems.  The  ultraviolet  (UV)  region  ranges  in  wavelength  from  0.2  to 
0.4  f an  .  The  visible  region  ranges  in  wavelength  from  0.4  to  0.7  pan.  The  near  infrared 
spectral  region  (NIR)  ranges  in  wavelength  approximately  from  0.7  to  1.1  pan  .  The  first 
infrared  band  is  the  short  wavelength  infrared  band  (SWIR)  which  approximately  covers 
1.1  to  2.5 pan.  The  second  infrared  band  is  the  mid- wavelength  infrared  (MWIR)  band 
and  covers  approximately  from  2.5  to  7.0  pan.  The  third  infrared  band  is  the  long 
wavelength  infrared  (LWIR)  spectral  band  and  it  covers  approximately  from  7.0  to 
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1 5  /jm.  The  fourth  infrared  band  is  the  far  infrared  (FIR)  band  or  very  long  wave  infrared 
(VLWIR)  band  which  covers  the  spectral  region  above  15  fjm.  [Ref.  5]. 

C.  REFRACTIVE  INDEX 

Due  to  interaction  that  occurs  between  the  wave  electric  field  and  the  medium  in 
which  it  propagates,  the  refractive  index  is  different  from  that  of  the  free  space.  We  can 
picture  the  electron  as  attached  to  the  end  of  a  spring  as  a  way  to  describe  this  interaction 
with  a  simple  electron  oscillator  model.  In  that  case  we  can  show  that  the  refractive  index 
is  a  complex  quantity  with  a  real  part  given  by: 

-  >  +  t  (8) 

j  \o)j -a>2)  +/jG)~ 

and  the  absorption  coefficient  (imaginary  part)  is  given  by: 

,  ,  w 

m£oc  j  [a2 -<v2)  +y-(Q- 

where  N  is  the  number  of  molecules  per  unit  volume 
q  is  the  charge  of  the  electron 
m  is  the  electron  mass 

co  is  the  frequency  that  drives  the  damped  harmonic  oscillator 
C0j  is  the  natural  radial  radian  frequency  of  the  electrons 

f  j  is  the  number  of  electron  with  frequency  co]  (oscillator  strength) 

Yj  is  the  damping  coefficient  in  each  molecule 
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The  real  part  of  the  refractive  index  represents  the  wave  dispersion  (i.e.,  controls  the 
phase  velocity)  while  the  imaginary  part  represents  wave  absorption  at  the  resonant 
frequency.  Close  to  a  resonance,  the  index  of  refraction  drops  sharply.  This  phenomenon 
is  called  anomalous  dispersion  because  of  this  atypical  behavior.  The  maximum 
absorption  occurs  in  that  region  and  that  happens  because  we  drive  the  electrons  to  their 
resonant  frequency.  The  result  is  that  a  very  large  amount  of  energy  is  dissipated.  Away 
from  the  resonances  the  index  of  refraction  simplifies  to  the  following  form: 


n-1 +^~  y,  /> 
2ms0  jjcoj -co  ) 


(10) 


Moreover,  for  transparent  materials,  the  nearest  significant  resonances  generally  lie  in  the 
ultraviolet  [Ref.  6],  so  that  co<  co] . .  Using  Taylor  series  we  can  show  that  in  that  case  Eq. 

10  becomes: 


n=l+ 


y2me0  j  co) j 


+  co 


Nq2  ^  fj 


2msoy  co  j 


(11) 


2  tic 

With  X  -  - we  can  write  Eq.  1 1  as 

co 


with 


A  =  Nl> 


2ms0-JcoJ 


and 


n=l+A(i+A) 

A 


(12) 


B  = 


V  03  j  J 
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This  is  Cauchy’s  equation  in  the  two  constants  form.  The  constant  A  is  called  the 
coefficient  of  refraction  and  the  constant  B  is  called  the  coefficient  of  dispersion. 

A  more  general  equation,  which  gives  the  index  of  refraction  over  a  large  wavelength 
range,  is  the  Sellmeier  equation,  which  is  also  valid  in  the  regions  close  to  resonances. 
This  equation  for  a  single  resonance  wavelength  A0  is  given  by  [Ref.  1] 


n2=l+ 


A  A2 
A2-!2 


In  cases  of  more  than  one  resonance  Eq.  1 3  becomes 


(13) 


2  ,  V-  A, A2 
n-l+ >  - 


(14) 


D.  REFRACTIVITY 

The  value  of  the  refractive  index  of  air  is  veiy  close  to  unity.  Sometimes  it  is 
convenient  to  discuss  refractive  effects  in  terms  of  refractivity.  The  refractivity  is 
proportional  to  the  difference  between  n  and  1  and  is  defined  as  R=(n-1)  xlO6  [Ref.  1, 
Ref.  7],  The  refractivity  of  the  atmosphere  through  the  visible  and  IR  is  shown  in  the 
Figure  5  [Ref.  1].  From  Figure  5  we  observe  that  in  the  visible  range  the  variation  with 
wavelength  is  small  and  in  the  near  infrared  is  negligible. 
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Refractivity  vs  wavelength 
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Figure  5.  Atmospheric  Refractivity  versus  wavelength.  Reported  measured  values  are 
compared  with  computation  using  the  Edlen  equation  (Eq.  15)  [Ref.l]. 


E.  ATMOSPHERIC  REFRACTION 

Due  to  the  electromagnetic  nature  of  many  sensors  used  by  today’s  military,  the 
EM  propagation  through  the  atmosphere  is  a  question  of  extreme  importance  and 
relevance.  The  atmospheric  refraction  arises  from  the  density  variation  of  air  which 
decreases  almost  exponentially  with  altitude.  So,  we  expect  that  fluctuations  in 
temperature  T,  humidity  Q  and  pressure  P  will  cause  fluctuations  in  both  the  real  and 
imaginary  parts  of  the  refractive  index.  Changes  in  the  vertical  gradient  and  to  a  lesser 
extend  horizontal  gradient  of  the  index  of  refraction  cause  the  atmospheric  refraction.  In 
wavelength  regions  remote  from  major  absorptions  the  refractive  index  can  be 
represented  for  dry  air  in  the  form  [Ref.  1]: 
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(15) 


(16) 


Where  A  is  the  wavelength  in  micrometers 
P  is  the  pressure  in  millibars  and 
T  is  the  temperature  in  degrees  Kelvin. 

In  the  IR  and  the  visible,  the  index  gradient  of  air  is  dominantly  dependent  on  the 
temperature  gradient  in  the  operational  windows  of  EO  systems,  where  absorption  is  low. 


F.  A  NON-REFRACTIVE  ATMOSPHERE 

In  a  non-refractive  atmosphere  (i.e.,  one  in  which  the  refractive  index  is  invariant 
with  position)  electro-optical  (EO)  radiation  follows  straight  line  trajectories  and  hence 
the  maximum  detection  ranges  are  limited  by  the  terrestrial  horizon.  The  range  limitation 
due  to  the  earth’s  curvature  is  called  horizon-limited  range  (HLR)  [Ref.  8]  and  for  a  given 
sensor  and  target  height,  can  be  obtained  using 

HLR  =V2a  (4hs+4ht)  x  10'3  s  3.57  (-JK  +  4K)  (17) 

where  HLR  is  expressed  in  kilometers 
hj  is  the  sensor  height  in  meters 

h,  is  the  target  height  in  meters  and 

a  is  the  mean  earth  radius  expressed  in  meters  (a  =  6370000  m) 
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This  value  for  the  mean  earth  radius  is  given  by  [Ref.  9,  Ref.  10].  Moreover  in  [Ref.  9] 
for  the  midlatitude  winter  atmosphere  and  custom  atmosphere'  the  radius  of  earth  (km) 
that  is  used  in  the  calculations  of  MODTRAN  is  a  =6371  km.  This  is  the  radius  of  earth 
at  the  particular  latitude  at  which  the  calculation  is  to  be  performed.  As  we  will  see  on 
Data  Analysis  the  atmospheres  that  will  be  used  for  the  IRTOOL  calculations  are  the 
midlatitude  winter  and  custom  atmosphere.  So,  Eq.  17  with  the  radius  of  earth  as 
specified  above  will  help  us  to  evaluate  the  non-refractive  calculations  of  IRTOOL. 


Figure  6.  A  non-refractive  and  a  refractive  atmosphere  with  the  respective  elevation 
angles  s°  and  e .  The  tangent  represents  the  observed  source  detection  [After  Ref  14]. 

G.  REFRACTION  EFFECTS 

The  surface  layer  of  the  marine  boundary  layer  influences  the  IR  sensing  of  low 
altitude  targets.  Three  atmospheric  effects  are  primarily  responsible  for  limiting  sensor 
performance  [Ref.  11]: 
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(a)  Ray  refraction 

(b)  Turbulent  distortion 

(c)  Atmospheric  extinction  by  molecules  and  aerosols 

In  this  chapter  we  will  describe  only  the  effects  of  the  ray  refraction  due  to  the 
vertical  variation  of  the  index  of  refraction  which  dominates  in  most  cases  the  horizontal 
variations.  The  sign  of  the  vertical  gradient  of  the  refractive  index  divides  the  refraction 
conditions  into  two  categories.  In  the  first  one,  the  gradient  is  positive  and  we  have 
subrefraction  (inferior  mirage)  and  in  the  second  one  the  gradient  is  negative  and  we  have 
superrefraction  conditions  (superior  mirage). 

Mirage,  as  defined  by  the  dictionary,  is  an  optical  illusion  caused  by  the 
atmospheric  conditions,  in  which  an  observer  sees  a  nonexistent  body  of  water  or  an 
image  of  some  object.  Examples  of  mirages  are  false  pools  of  water  seen  over  hot  desert 
sands  or  hot  pavement,  and  at  sea  an  inverted  image  of  a  ship  seen  in  the  sky.  These 
phenomena  are  well  understood  and  can  be  explained  by  two  facts.  First  the  light  rays  are 
bent  when  they  pass  between  media  of  different  densities  and  second,  the  boundary 
between  two  such  media  acts  as  a  mirror  for  rays  of  light  coming  in  at  certain  angles  [Ref. 
12].  Mirages  can  be  photographed  (Figure  7).  For  most  mirage  effects  variations  in  the 
temperature  mainly,  but  also  in  the  pressure  and  the  water  vapor  content  of  the 
atmosphere,  cause  corresponding  variation  in  the  refractive  index  that  causes  the  bending 
of  light  rays  as  is  predicted  by  Snell’s  law  of  refraction.  The  optical  ray  bending  is 
determined  primarily  by  the  temperature  gradient  [Ref.  13] 
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dT  _  Z 
dz  akz 


(18) 


where  T*  is  a  temperature  scaling  parameter,  a  is  the  ratio  of  the  heat  transfer  to 
momentum  transfer  and  k  is  the  von  Karman  constant,  approximately  equal  to  0.4.  At 

low  elevation,  the  term  -  in  Eq.  18  results  in  a  very  large  temperature  gradient  and 

Z 

excessive  ray  bending.  Before  describing  in  more  detail  the  subrefraction  and 
superrefraction  conditions  we  should  have  in  mind  that  a  light  ray  always  bends  toward 
the  denser  medium  (colder  air). 


Figure  7.  The  desert  (or  hot)  road  mirage  [From  Ref.  3]. 
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H.  SUBREFRACTION  (INFERIOR  MIRAGE) 


The  kind  of  mirage  formed  over  a  hot  surface  is  sometimes  referred  to  as  an 
inferior  mirage. 


Figure  8.  Ray  diagram  for  subrefraction  conditions  [From  Ref.  8]. 


The  term  inferior  mirage  refers  to  the  fact  that  the  upright  mirage  appears  below  its  real 
position.  This  effect  occurs  under  conditions  of  negative  Air  to  Sea  Temperature 
Difference  (ASTD).  From  Figure  8  we  can  see  how  the  light  rays  propagate  through  a 
subreffactive  atmosphere.  Rays  leaving  the  source  at  angles  below  the  horizon  will  be 
refracted  upwards  due  to  the  refractive  gradient  leading  to  a  crossing  of  ray  trajectories. 

The  following  three  Figures  9,  10,  11  will  help  us  to  understand  the  results  given 
from  the  data  analysis  in  Chapter  VII. 
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Figure  9.  Ray  paths  under  subreff  active  conditions  [From  Ref.  3]. 

In  the  desert  the  sand  becomes  hotter  than  the  air  and  heats  the  air  layer  immediately 
above  the  surface  to  a  higher  temperature  than  the  air  a  few  meters  above.  The  greatest 
bending  of  the  ray  occurs  closest  to  the  ground  because  the  temperature  there  is  changing 
more  rapidly  with  the  height.  In  Figure  9  we  see  five  rays  coming  from  the  back  of  the 
monkey  part  way  up  the  tree  trunk.  The  rays  start  off  in  different  directions  and  hit  a  wall 
at  a  specific  range  from  the  tree.  As  we  can  see,  rays  4  and  5,  although  initially  directed 
below  the  other  three,  end  up  higher  on  the  wall.  If  we  put  our  eye  at  A  we  will  see  two 
views  of  the  monkey,  one  direct  and  the  other  inverted.  The  same  happens  at  point  B  only 
now  the  two  images  are  closer  than  before.  At  point  C  the  two  images  have  become  one, 
and  below  C  we  do  not  see  any  image.  The  point  C  is  the  lowest  point  at  which  a  ray  can 
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be  received  at  this  particular  distance,  as  shown  in  Figure  10.  With  our  eye  at  point  C,  we 
see  the  top  of  the  tree  but  none  of  the  tree  trunk  below  the  monkey. 


Figure  10.  The  origin  of  a  vanishing  line  in  the  desert  mirage  [From  Ref.  3]. 


Figure  11.  "The  appearance  of  the  desert  mirage  scene  with  (A)  the  observer  positioned 
as  shown  in  Figure  above  and  (B)  the  observer  further  from  the  tree"  [From  Ref.  3]. 
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Figure  1 1  (A)  shows  the  image  we  see  looking  from  point  C.  There  is  a  vanishing  line 
passing  through  the  monkey’s  back.  Every  point  on  the  tree  below  that  line  cannot  be 
seen  but  every  point  above  that  line  has  two  images  and  one  of  those  two  images  is 
inverted.  If  we  move  further  away  from  the  tree  the  vanishing  line  will  rise  as  shown  in 
case  B.  That  happens  because  all  the  rays  will  pass  above  our  head.  Finally,  if  we  move 
even  further  the  tree  will  disappear.  From  Figures  8,  12  we  see  that  subrefraction 
conditions  impose  an  absolute  limit  of  detection  range  shorter  than  the  terrestrial  horizon 
(HLR). 


Figure  12.  Under  subrefractive  conditions  two  images  of  the  source  are  seen  at  two 
elevation  angles.  L°h  and  s°h  are  the  horizon  distance  and  angle  [From  Ref.  14]. 
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In  coastal  environments  the  appearance  of  subrefraction  conditions  is  more  common, 
especially  when  those  conditions  are  driven  by  the  sea  wind.  For  those  conditions  the 
expression  of  the  maximum  intervision  range  is  used  to  “denote  the  absolute  detection 
range  limit  imposed  by  refraction  for  a  given  sensor  and  target  height”  [Ref.  8].  This 
range  corresponds  to  the  HLR  for  non-refractive  conditions. 

I.  SUPERREFRACTION  (SUPERIOR  MIRAGE  OR  DUCTING  OR 
LOOMING) 

This  effect  arises  when  warmer  air  moves  over  a  colder  lake  or  ocean.  That  means 
that  we  have  positive  ASTD  conditions  causing  light  rays  to  be  bent  toward  the  earth  and 
hence  making  an  object  appear  to  be  raised  above  its  true  position  (the  object  is  seen 
above  its  actual  position).  When  this  condition  is  well  developed  we  can  see  objects  such 
as  a  distant  ship  beyond  the  horizon  or  the  sun  after  it  has  set  (Figure  13).  For  a  long  time 
this  phenomenon  has  been  called  looming  by  sailors.  The  superrefractive  conditions  are 
shown  in  Figures  14,  15.  We  see  that  under  these  conditions,  the  downward  bending  of 
the  rays  can  cause  the  detection  of  targets  well  beyond  the  horizon  even  though  the  target 
is  not  observable  under  non-refractive  conditions  (straight  line  propagation). 

These  conditions  are  not  so  often  encountered.  Depending  on  the  season  and  the 
region  of  the  world  it  is  possible  for  them  to  occur  near  the  coast,  especially  when  warm 
continental  air  from  the  land  is  blown  over  cool  water. 
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Figure  13.  Rays  producing  a  superior  mirage  [From  Ref.  3].  The  observed  object  is  seen 
as  an  elevated,  erect  image.  The  object  may  be  over  the  geometrical  or  "normal 
refraction"  horizon. 


Figure  14.  Under  superrefractive  conditions  a  target  can  be  seen  even  when  it  is  below 
the  horizon.  is  the  elevation  angle  [Ref.  14], 
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Figure  15.  Ray  diagram  for  superrefraction  conditions  [From  Ref.  8].  As  we  can  see  the 
downward  bending  of  the  rays  can  cause  the  detection  of  targets  well  beyond  the  horizon 
even  though  the  target  is  not  observable  under  non-refractive  conditions.  The  dot  line 
crossing  the  ray  trajectories  corresponds  to  HLR. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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III.  IRBLEM 


A.  GENERAL  DESCRIPTION 

The  IRBLEM  (IR  Boundary  Layer  Effects  Model)  is  a  software  package, 
consisting  of  various  modules,  which  can  be  used  as  executable  stand-alone  programs.  It 
is  developed  by  the  Defense  Research  Establishment,  Valcartier,  (DREV)  Canada,  and 
performs  calculation  of  optical  properties  in  a  propagation  volume  defined  by  the  user. 
Atmospheric  effects  such  as  refraction  and  turbulence  in  the  marine  surface  layer  are  also 
considered  in  contrast  with  MODTRAN. 

The  surface  layer  theory  used  by  IRBLEM  is  in  general  valid  for  the  first  30 
meters  above  the  surface  but  with  direct  calls  to  MODTRAN  calculations  can  be  made 
for  heights  above  this  defined  surface  layer. 

The  IRBLEMv3.1  code  was  downloaded  from  the  DREV  ftp  site  and  is  a  PC 
executable  program  compatible  only  with  the  MODTRAN  3vl.5.  Figure  16  shows  the 
internal  modular  structure  of  IRBLEM  with  the  meteorological  inputs  situated  at  the  top 
and  the  various  outputs  at  the  bottom.  As  we  can  see,  in  a  surface  layer  domain  defined 
by  the  user,  vertical  profiles  of  refractivity  M(h),  of  aerosol  extinction  coefficients  a(h), 
and  of  the  refractive  index  structure  parameter  C 2  (h)  are  calculated  in  the  first  phase. 

Also,  depending  on  whether  or  not  MODTRAN  is  used,  the  molecular  transmittance 
versus  height  is  available.  In  the  second  phase,  all  the  above  mentioned  profiles  are  used 
as  inputs  in  the  ray  tracing  program,  which  calculates  the  optical  properties  and 
atmospheric  effects. 
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Figure  16.  IRBLEM  functional  diagram  [From  Ref.  15]. 


B.  MAIN  INPUTS 

The  main  inputs  can  be  divided  into  five  groups.  The  first  group  is  the 
propagation  domain,  the  second  the  meteorological  conditions,  the  third  the  sensor  data, 
the  fourth  the  running  options  and  the  fifth  one  the  ray  tracing  parameters.  The  active 
parameters  that  are  used  for  each  group  are  shown  in  the  following  Tables  1-5. 
PROPAGATION  DOMAIN  DATA  (GROUP  1) 

Minimum  Range  and  Maximum  Range  (km) 

Maximum  height  (m) 

Range  sampling/step  (km) 

Height  sampling/step  (m) 

Table  1.  Propagation  domain  parameters  in  IRBLEM  [After  Ref.  15]. 
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METEOROLOGICAL  DATA  (GROUP  2) 

Height  of  the  Air  Temperature  sensor,  Pressure  sensor,  and  Humidity  sensor  in  meters 

Wind  speed  (m/sec)  and  height  of  the  wind  speed  measurement  in  meters 

24  hr  average  wind  speed  (m/sec) 

Pressure  (mbar) 

Air  Temperature  (deg  C) 

Water  Temperature  (deg  C) 

Humidity  (%) 

Wave  height  (H  X )  (m) 

Visibility  (Km) 

Solar  Radiation  (W/m 2 ) 

Air  Mass  parameter  with  values  between  1-10 

Precipitation  type  (rain,  drizzle,  snow,  fog,  no  weather) 

Precipitation  rate  (mm/h) 

Fog  liquid  water  content  (g/m3) 

Table  2.  Meteorological  data  active  inputs  in  IRBLEM  [After  Ref.  15]. 

SENSOR  DATA  (GROUP  3) 

Waveband  (min  and  max  wavelength)  in  ( /u  m) 

Sensor  Height  (m) 

Table  3.  Active  sensor  parameters  that  are  used  by  IRBLEM  [After  Ref.  15]. 
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RUN  OPTIONS  (GROUP  4) 

Job  identification  number  (1-999) 

Spectral  resolution  (cm  ) 

Aerosol  computation  mode  (NOV AM  or  Kel) 

Use  of  MODTRAN  (yes  or  no) 

Meteorological  profiles  mode  (“wavy”  or  “non-wavy”) 

Table  4.  Running  options  parameters  in  IRBLEM  [After  Ref.  15]. 

RAY-TRACING  PARAMETERS  (GROUP  5) 

Minimum  ray  launching  angle  (deg) 

Maximum  ray  launching  angle  (deg) 

Maximum  number  of  rays  that  will  be  used  (integer) 

Nominal  ray  angular  separation  (deg) 

Nominal  ray  integration  step  (m) 

Minimum  ray  separation  angle  (deg)-when  approaching  horizon 

Range  integration  step  for  path  radiance  calculation  (km) 

Wave  cut-off  factor  (unitless) 

Total  extinction  surface  sublayer  (m) 

Table  5.  Ray  tracing  parameters  in  IRBLEM  [After  Ref.  15]. 


For  the  evaluation  of  the  sea-particle  concentration  as  function  of  the  wind  speed 
two  models  are  available: 
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(1)  The  NOV AM  relationships 

(2)  A  new  model  developed  by  Kel  Research 

Also,  for  the  computation  of  the  reffactivity  profile  we  have  two  options:  the  “wavy”  and 
the  “non-wavy”.  In  the  first  one  the  sea  wave  effect  (the  wave  structure)  is  taken  into 
account  and  modifies  the  reffactivity  profile  produced  by  the  non-wavy  mode.  The  profile 
produced  by  the  non-wavy  mode  is  averaged  out  over  a  sea-wave  period  as  a  function  of 
the  wave  height.  According  to  Dion  [Ref.  15]  “it  is  deemed  preferable  to  use  the  non- 
wavy  mode  until  validation  has  been  achieved.” 

Something  that  will  play  a  crucial  role  to  the  data  analysis  is  the  FWH  (Factor 
Wave  Height).  The  ray  cut-off  height  due  to  waves  is  given  in  IRBLEM  by  the  following 
expression: 

Cut-off  height  =  (FWH  xH,)  +  minabslayer  ( 1 9) 

That  means  from  the  water  surface  to  the  cut-off  height  the  waves  are  assumed  to  produce 
a  fully  absorbing  layer  that  stops  all  the  rays.  The  H  w  or  H  is  the  significant  wave 

height  in  meters  and  minabslayer  is  a  sublayer  above  the  surface  where  total  extinction  is 
assumed  to  exist.  The  significant  waveheight  is  taken  here  to  equal  four  times  the  rms 
waveheight  [Ref.  11,  Ref.  13].  According  to  [Ref.  15]  for  the  ray-tracing  parameters 
These  parameters  are  not  likely  to  require  modifications  on  a  run-to-run  basis  under 
normal  use  of  IRBLEM.  Default  values  should  be  adequate  under  most  circumstances.” 
The  default  value  for  the  FWH  factor  is  given  as  0.6.  So,  we  expect  that  this  factor  will 
produce  a  significant  wave  height  much  smaller  than  the  measured  one  (about  50% 
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smaller).  In  that  case  we  expect  the  beginning  of  the  mirage  range  to  be  at  shorter 
distances. 


C.  MODULES  AND  UTILITY  PROGRAMS 

The  main  programs  and  utility  programs  that  compose  the  IRBLEM  v3.1  package 

are  shown  in  Tables  6,  7  [Ref.  15].  The  main  programs  with  the  asterisk(*)  are  the 
programs  that  are  called  by  IRBLEM. 


PROGRAM 

PURPOSE 

IRBLEM 

Main  program  —  Batch  or  Menu  —  based  interactive  mode 

LWWKD* 

Computation  of  modified  refractivity  profile,  M(h) 

WKDCN2* 

Computation  of  C  \  (h),  the  index  structure  parameter  profile 

WKDAERX* 

Computation  of  spectral  aerosol  extinction  coefficient  versus  height, 
aa(A,h) 

TRANSM* 

Calculation  of  molecular  transmittance  (using  MODTRAN  when 
available) 

RTR* 

Ray  -  tracing  program 

IRBLEMPP 

IRBLEM  post  processor 

BKGD 

Give  background  radiance  in  the  marine  surface  layer  _ 

Table  6.  Main  programs  of  IRBLEM  [From  Ref.  15]. 


UTILITY 

PURPOSE 

MODSET 

Set  path  to  MODTRAN  (when  available) 

MERGPROF 

Determine  atmospheric  profile  for  calculation  above  the  surface  layer 

MSGCHECK 

To  check  for  warnings  and  error  messages 

EXTR 

Extraction  of  results  (scalar  vs.  height,  spectral  vs.  wavelength)  for 
display 

CLEANER 

Delete  inputs/outputs  for  a  given  IRBLEM  job  ID 

Table  7.  Utility  programs  of  IRBLEM  [From  Ref  15]. 
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Now,  lets  describe  in  a  more  little  details  those  modules  that  were  used  in  this 

thesis. 

1.  LWWKD 

The  LWWKDv6.5  is  a  bulk  aerodynamic  model  developed  by  DREV  that 
produces  the  vertical  refractivity  profile  and  applies  to  RF,  IR  and  visible  wavelengths.  It 
is  based  on  the  Monin-Obhukov  similarity  theory.  This  model  uses  standard 
meteorological  measurements  [Ref.  15]  as  inputs,  namely: 

Water  temperature  (deg  C) 

Air  temperature  (deg  C) 

Air  humidity  (%) 

Wind  speed  (m/sec) 

Pressure  (mbar) 

Also  as  input  is  used  the  wavelength  and  the  choice  of  “wavy”  (LWWKD)  or  “non-wavy” 
mode  (LWKD).  Actually,  this  module  produces  the  vertical  profile  of  modified 
refractivity  up  to  50  meters  above  the  water.  The  modified  refractivity  is  defined  as  [Ref. 
15,  Ref.  23] 

M{h)  =  N(h)  +  0A57h  (20) 

where  N  =  (n-l)xlO6  is  the  refractivity,  n  is  the  index  of  refraction,  h  is  the  height  in 
meters  and  0.157  is  a  constant  in  units/m.  In  addition  to  the  vertical  profile  of  modified 
refractivity,  LWWKD  produces  the  profiles  of  its  first  and  second  derivative  and  also  the 
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vertical  profiles  of  temperature,  relative  humidity,  pressure  and  wind  speed.  The 
expressions  for  the  first  and  second  derivative  of  the  modified  refractivity  are  given  by: 


-  dNW.  +  q  1 5  (first  derivative)  (2 1 ) 

dh  dh 


d^Mjh)  _  d^Njh) 
dh 2  dh2 


(second  derivative) 


(22) 


2.  RTR 

The  RTRv6.2  is  a  program  developed  by  DREV,  which  for  a  user-defined  marine 
surface  propagation  layer  integrates  the  atmospheric  properties.  In  that  calculation,  ray 
paths  calculated  from  the  sensor  (height=sensor  height  and  range=0  km)  depend  on  the 
vertical  profiles.  The  Earth’s  curvature  also  is  taken  into  account. 

The  RTR  uses  as  inputs  the  outputs  of  all  the  other  modules.  The  outputs  of  this 
program  are  many  interesting  quantities  such  as  the  MIVR,  the  geometrical  horizon,  the 
beginning  of  the  mirage  zone,  the  angular  deviation  i.e.,  the  angular  difference  between 
the  elevation  of  the  initial  ray  and  the  elevation  of  the  initial  ray  that  hits  the  same 
coordinates  under  non-refractive  conditions,  the  geometry  of  every  ray,  the  refractance, 
aerosol  transmittance,  spectral  path  radiance  and  transmittance  along  each  ray  and  many 
more. 

Here,  we  can  say  that  the  MIVR  and  also  the  beginning  of  the  mirage  zone  are 
meaningful  quantities  only  under  unstable  conditions  (subrefraction  conditions-air  cooler 
than  water).  In  that  case  for  a  constant  altitude  target,  a  mirage  appears  between  the 
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beginning  of  the  mirage  zone  range  and  the  MIVR.  Under  stable  conditions 
(superrefraction  conditions-air  warmer  than  water)  mirage  does  not  occur.  In  that  case  the 
concept  of  MIVR  does  not  apply,  because  the  rays  can  reach  very  long  ranges  near  the 
surface. 

3.  IRBLEMPP  (IRBLEM  POST  PROCESSOR) 

This  is  a  program  that  describes  the  cumulative  atmospheric  optical  effects  on 
imaging  of  a  point  target,  which  is  located  at  coordinates  given  by  the  user  (height, 
range). 

For  a  given  target  position  the  IRBLEMPP  calculates  [Ref.  15]  the  apparent 
elevation  of  horizon,  the  elevation  of  the  geometrical  horizon,  the  minimum  observable 
height  (the  minimum  height  above  the  water  where  a  target  can  be  detected  at  a  given 
distance),  the  apparent  target  position,  the  refractance,  information  on  the  mirage  image 
(if  a  mirage  is  present)  as  well  as  for  the  direct  (real)  image,  spectral  quantities  such  as 
molecular  and  aerosol  transmittance,  path  radiance  and  many  other  parameters. 
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IV.  IRTOOL 


A.  GENERAL  DESCRIPTION 

IRTOOL  is  an  IRST  X-windows  analysis  tool,  which  was  developed  by  Arete 
Associates  and  NSWC  under  the  sponsorship  of  the  Office  of  Naval  Research  in  support 
of  the  Infrared  Analysis  Modeling  and  Measurements  Program  (IRAMMP).  It  is  oriented 
towards  IRST  analysis  [Ref.  16].  This  computer  program  is  a  simulation/computational 
model  environment,  which  uses  a  number  of  stand-alone  executable  modules  in  a  unified 
package,  to  produce  both  a  visual  simulation  and  engineering  calculations  for  scenarios 
appropriate  to  the  Infrared  Search  and  Track  system. 

IRTOOL  includes  a  variety  of  sensor  performance  models,  a  few  point  source  IR 
signature  models,  an  atmospheric  propagation  model  (MODTRAN  and  LOWTRAN) 
including  standard  atmospheric  models  and  provision  to  build  a  local  atmospheric  model, 
a  turbulence  model  and  an  atmospheric  refraction  model.  The  entire  package  is  designed 
so  that  the  users  may  insert  local  modules  for  their  purposes.  IRTOOL  applies  to  the  low 
altitude  ocean  boundary  layer  regime.  Other  features  of  IRTOOL  other  than  its  ability  of 
producing  simulated  ocean  and/or  sky  scenes  including  clouds  is  the  injection  of  targets 
into  both  simulated  and  real  scenes.  This  level  of  analysis  allows  visualization  of  the 
fluctuations  in  the  observed  target  signal  due  to  the  obscuration  by  the  wave  crests. 

For  this  thesis  we  used  IRTOOL  version  2.1.0  which  is  installed  on  the  Sun 
machines  in  the  Simulation  Laboratory  of  the  Physics  Department  in  the  Naval 
Postgraduate  School  (NPS). 
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B.  INPUT  PARAMETERS 

The  IRTOOL  parameters  are  organized  by  subject  and  each  subject  includes  one 
or  more  group  of  parameters.  Table  8  lists  the  group  of  parameters  for  each  subject. 

The  group  of  the  input  parameters  is  found  under  the  Input  menu.  The  user  has 
three  options:  (1)  load  a  previously  saved  set  of  parameters,  (2)  edit  the  present  parameter 
values,  (3)  load  a  default  set  of  parameter  values  and  (4)  use  Advanced  flow  control. 
The  fourth  option  is  used  when  the  user  is  doing  several  similar  calculations  with  similar 
inputs.  In  this  case,  previously  calculated  results  for  the  different  runs  can  be  reused  and  a 
significant  amount  of  time  can  be  saved.  It  was  very  useful  in  this  thesis  to  use  that 
option,  especially  when  we  had  to  use  the  same  LOWTRAN  file  for  a  number  of  different 
runs.  The  most  important  consideration  in  using  the  advanced  flow  control  is  to  ensure 
that  the  inputs  are  consistent  with  the  results  we  are  using. 

C.  MODULES 

IRTOOL  consists  of  eighteen  modules  which  are  stand  alone  executable  programs 
(or  groups  of  programs)  with  well-defined  inputs,  outputs  and  algorithms.  The  sequence 
in  which  those  modules  are  executed  is  very  important  because  some  modules  require  as 
inputs  some  of  the  outputs  from  other  modules.  Module  communication  is  accomplished 
with  disk  files.  A  brief  description  of  the  modules,  the  parameter  group  that  is  used  for 
each  of  them  and  also  their  dependencies  are  shown  in  Tables  9,  10.  For  a  specific  output 
one  or  a  number  of  modules  must  be  executed. 
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Subject 

Parameter 

Group 

Comments 

Sensor 

Type 

Contains  the  sensor  system  parameters  such  as  the  detector 
IFOY,  detector  spacing  (also  known  as  pitch  or  sampling), 
the  sensor  passbands,  and  detector  efficiencies 

Polarizer 

Input  parameters  which  describe  the  polarizing  element,  if 
any,  attached  to  the  sensor. 

Operation 

Contains  the  input  parameters  which  specify  where  the 
sensor  is  looking,  at  what  altitude  the  sensor  is  located,  and 
how  the  sensor  platform  is  moving. 

Scene  effects 

Input  parameters  which  affect  how  the  simulated  sensor 
response  is  applied  to  an  image. 

Target 

Type 

Target  size  and  emissivity. 

Trajectory 

Input  parameters  which  control  the  position  and  velocity  of 
the  target  including  the  target  speed,  altitude,  and  ranges  of 
interest.  This  menu  has  5  different  forms  depending  on  how 
the  trajectory  is  specified.  The  Pursuit/Vectored  option  is  for 
targets  with  constant  speed  and  altitude  on  collision  courses 
with  the  sensor.  With  the  Specify  Range/Bearing  and 
Range/Bearing  from  File  options  the  user  defines  the  target 
position  relative  to  the  sensor  as  a  function  of  time.  With  the 
Specify  Speed/Heading  Option  and  Speed/Heading  from  File 
options  the  user  provides  an  initial  target  position  and  defines 
the  target  velocity  as  a  function  of  time.  The  last  four  options 
allow  variable  target  speeds  and  altitudes. 

Injection 

Input  parameters  which  affect  how  the  simulated  target  will 
be  injected  into  an  image. 

Environment 

Atmosphere 

Input  parameters  which  control  the  LOWTRAN 
atmosphere,  solar  position  and  how  MODTRAN 
calculations  are  performed. 

Ocean 

Input  parameters  describing  physical  properties  of  the  ocean 
such  as  swell  and  parameters  which  affect  how  the  ocean 
surface  is  modeled. 

Cloud 

Input  parameters  specifying  clear-sky  or  a  cloud  type  and  the 
associated  cloud  parameters 
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Scene 

Simulator 

Scene 

Simulator 

Input  parameters  which  specify  the  number  and  size  of  the 
images  to  be  generated  by  the  scene  simulator  module. 

System 

Performance 

Performance 

Input  parameters  which  set  the  number  of  analytic  system 
performance  calculations  to  perform  and  the  assumptions 
made  about  the  type  of  processing  to  model. 

Data 

Processing 

IRAMMP 

Format 

Converter 

Input  parameters  for  the  program  which  converts  image  data 
to  HDF  format. 

Stream 

Processing 

Input  parameters  which  control  the  input  and  behavior  of 
a  “stream”  of  modules,  created  in  a  run  script. 

Preprocessor 

Input  parameters  for  controlling  image  processing  module. 
Processing  options  include  bias  removal,  detrend, 
subpatching. 

Statistical 

Analysis 

Input  parameters  for  the  program  which  does  statistical 
analysis  (mean,  variance,  pdf,...)  on  images 

Spectral 

Analysis 

Input  parameters  for  the  program  which  does  spectral 
analysis  (power  spectral  density,  coherence)  on  images. 

User 

User 

Parameters  for  user  supplied  modules. 

Several  sample  user  modules  are  included  with  IRTool,  and 
their  input  parameter  sets  are  available  to  load. 

Advanced 

Flow  Control 

Advanced  Flow 
Control 

Tells  modules  where  to  find  inputs.  Use  reserved  for 
experienced  users  interested  in  reusing  previously  calculated 
results  in  order  to  save  on  run  time. 

Table  8.  Input  parameters  by  category  in  IRTOOL  [After  Ref.  16]. 


The  user  has  the  capability  to  see  what  modules  must  be  run  and  in  what  order  for  a 
specific  output,  by  selecting  between  two  options.  Those  are  the  Default  Run  Script  and 
the  Advanced  Run  Script  under  the  Run  option  menu.  These  options  and  also  for  many 
other  capabilities  of  IRTOOL  will  be  referred  to  in  the  following  paragraphs. 
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In  this  thesis  we  used  some  of  the  modules  but  not  all  of  them.  We  will  describe 
in  a  little  more  detail  the  modules  that  were  used,  the  specific  inputs  for  those  modules 
and  also  their  outputs. 

1.  Atmosphere  Profiler  Module 

This  module  generates  the  vertical  profiles  of  temperature,  humidity,  wind  speed, 
refractive  index  structure  parameter,  pressure,  mixing  ratio  and  many  other  parameters. 
The  calculation  of  those  profiles  is  performed  with  the  use  of  the  bulk  aerodynamic 
method  and  the  similarity  theory  of  Monin-Obkhov  [Ref.  16].  In  order  to  execute  this 
module  we  need  to  enter  a  specific  group  of  input  parameters.  Those  parameters 
correspond  to  the  group  parameter  Environment,  and  are  shown  in  Table  1 1 . 

As  we  can  see  from  the  input  parameters  in  Table  11,  there  is  no  input  for 
pressure  measurement.  IRTOOL  also  assumes  that  all  the  measurements  of  Temperature, 
Humidity  and  Wind  speed  are  made  at  the  same  height.  This  is  a  deficiency  of  the 
program  in  contrast  with  IRBLEM,  in  which  for  each  measurement  we  must  specify  the 
height  at  which  the  measurement  was  made. 

When  IRTOOL  performs  MODTRAN  calculations  it  does  two  things.  First,  it 
generates  a  MODTRAN  input  file  made  either  from  the  input  menus  or  a  combination  of 
the  menus  with  a  custom  input  file  provided  from  the  user.  Then  using  this  profile,  it 
performs  a  sequence  of  MODTRAN  calculations  and  produces  the  vertical  profiles  of  the 
atmospheric  parameters.  We  say  vertical  because  is  assumed  that  we  do  not  have 
horizontal  variations.  Also  from  Table  11  we  see  that  we  have  the  option  to  select  a 
custom  atmosphere.  In  this  case  the  user  must  provide  his  own  LOWTRAN  file  with 
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measurements  of  temperature,  pressure,  and  humidity  versus  height.  Those  measurements 
will  transform  into  a  LOWTRAN  file  that  will  be  the  basis  for  the  generation  of  the 
vertical  profiles  affecting  the  refractivity. 


Module 

Comments 

Sensor  Model 

Calculation  of  sensor  transfer  function  and 
noise 

Engagement 

Position  of  target  relative  to  sensor  without  ray 
bending 

Atmosphere  profiler 

Profile  of  atmospheric  boundary  layer  and 
Generation  of  its  own  LOWTRAN  file 

Atmosphere  without  radiance 

Calculation  of  refraction  (ray  paths  to 
background  and  target)  and  atmospheric 
distortions  for  the  scene  and  for  the  target,  if 
present 

Atmosphere  with  radiance 

Same  as  Atmosphere  without  radiance  but  also 
calculates  path  radiances  and  transmissivities 
to  background  and  target 

Target  Scintillation 

Detailed  calculation  of  target  scintillation  due 
to  atmospheric  turbulence 

Sky  Radiance 

Calculates  a  sky  radiance  dome  as  viewed  from 
the  surface  of  the  ocean 

Cloud  Clutter  for  ESNR 

Calculation  of  cloud  PSDs  at  points  in  target 
trajectory 

Scene  Simulator 

Generate  images  of  background  and 

Elevation  profiles  of  the  mean  and  variance 

Target  Signature 

Signature  image  at  target  range  and  look  aspect 
angle 

Target  Injector 

Output  scenes  with  injected  target  signature 

Sensor  Effects 

Add  effects  of  noise,  jitter,  and  system  transfer 
function  to  simulated  images 

System  Performance 

Analytic  estimate  of  signal  to  noise  ratio 

IRAMW  Data  Format  Conversion 

Conversion  of  NATO  format  to  HD 

Preprocessor 

Variety  of  options  for  orienting,  subpatching, 
and  processing  HDF  data 

Statistical  Analysis 

Mean,  variance,  correlation  and  probability 
density  estimates  from  HDF  data 

Spectral  Analysis 

PSD  and  cross  spectral  density  estimate  from 
HDF  data 

Table  9.  IRTOOL  module  description  [After  Ref.  16]. 
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Module 

Parameter  Group  Used 

Other  Modules  to  Run 

Sensor  Model 

Sensor  Type 

Sensor  Polarizer 

None 

Engagement 

Sensor  Operation 

Target  Trajectory 
(Sensor  Type  and 

Scene  Simulation  used 
for  diagnostics) 

None 

Atmosphere  Profiler 

Environment 

None 

Atmosphere  Effects 

Sensor  Type 

Sensor  Operation 

Environment 

Atmosphere  Profiler 
Engagement 

Cloud  Clutter  for  ESNR 

Sensor  Type 

Sensor  Operation 

Environment 

Target  Trajectory 

Atmosphere  Profiler 
Engagement 

Atmosphere  Effects 

Sky  Radiance 

Sensor  Type 

Sensor  Operation 

Environment 

Atmosphere  Profiler 
Atmosphere  Effects 

Scene  Simulator 

Sensor  Type 

Sensor  Operation 

Environment 

Scene  Simulation 

Atmosphere  Profiler 
Atmosphere  Effects 

Sky  Radiance  (if  ocean  is  in 
FOV) 

Target  Signature 

Sensor  Type 

Atmosphere 

Target  Type 

Target  Trajectory 

Engagement 

Atmosphere  Profiler 
Atmosphere  Effects 

Scintillation 

Sensor  Type 

Sensor  Operation 

Target 

Atmosphere 

Atmosphere  Profiler 
Atmosphere  Effects 

Engagement 

Target  Injector 

Sensor  Type 

Sensor  Polarizer 

Sensor  Operation 

Atmosphere 

Sensor  Model 

Engagement 

Atmosphere  Effects 

Scintillation 

Scene  Simulator 

Data  Converter 

Data  Converter 

None 

Sensor  Effects 

Sensor  Type 

Sensor  Polarizer 

Sensor  Scene  Effects 

Sensor  Model 

Scene  Simulator 

Target  Injector 

Table  10.  Inputs  used  by  each  module  by  Parameter  Group  and  Module  [After  Ref.  16]. 
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PARAMETERTS 

COMMENTS 

LOWTRAN  Atmosphere  Type 

The  user  can  select  one  of  the  six  available 
model  atmospheres  for  the  generation  of  the 
LOWTRAN  file.  These  Atmospheres  are  the 
U.S.  Standard  Atmosphere,  Tropical 
Atmosphere,  Midlatitude  Summer  Atmosphere, 
Midlatitude  winter  Atmosphere,  Subarctic 
Summer  Atmosphere,  Subarctic  Winter 
Atmosphere  and  the  Custom  Atmosphere. 

Solar  Position 

Options:  Night  or  Day.  If  we  select  Day  then 
we  have  to  specify  the  Solar  Azimuth  and  also 
the  Solar  Elevation 

Constant  Flux  Surface  Layer 

Options:  On  or  Off.  If  we  choose  On  then  we 
produce  a  constant  flux  layer  with  depth 
specified  as  Surface  Layer  Height  under  the 
Algorithm  Control  dialog  box 

Reference  Height  in  meters 

The  height  above  the  mean  sea  level  in  which 
measurements  of  temperature,  relative 
humidity  and  wind  speed  were  made.  Must  be 
between  1-30  meters. 

Air  Temperature  at  the  reference  height 

In  degrees  Celsius 

Relative  Humidity  at  the  reference  height 

In  percentage  (%) 

Wind  speed  at  the  reference  height  and  24hr 
average  wind  speed 

In  m/sec 

Wind  heading 

In  degrees 

ASTD  or  Water  Temperature 

In  degrees  Celsius 

Algorithm  Control 

We  specify  the  surface  layer  height  refraction, 
the  LOWTRAN  spectral  resolution  and  also 
we  select  the  coastal  influence  parameter 
(ICSTL)  between  the  values  1-10  and  whether 
we  will  use  refraction  or  straight  geometry. 

Ocean  parameters 

Swell  parameters  such  as  rms  amplitude,  peak 
frequency,  heading  etc. 

Cloud  parameters 

Options:  Clear  Sky  or  Cloud  type.  If  we  select 
the  second  one  we  have  to  enter  the  type  of  the 
clouds,  the  cloud  base  and  top  altitude,  the 
heading  and  speed. 

Table  1 1 .  Required  input  parameters  for  the  Atmosphere  Profiler  Module  [After  Ref.  16]. 
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In  this  thesis  we  used  for  each  data  set  two  atmospheres:  The  midlatitude  winter 
atmosphere  and  the  custom  atmosphere,  because  we  would  like  to  investigate  how  the 
absence  of  the  pressure  entry  in  input  parameters  affects  the  results.  A  second  reason  was 
to  see  if  the  "deficiency"  of  IRTOOL  to  use  only  one  reference  height  for  all  the 
measurements  can  alter  significantly  the  reffactivity  calculated  from  IRTOOL  results.  For 
those  reasons  we  ran  for  the  same  entries  the  NPS  Profiler  (a  bulk  aerodynamic  model  to 
which  we  will  refer  later)  and  we  used  its  predicted  vertical  profiles  of  temperature, 
pressure  and  relative  humidity  as  inputs  for  the  generation  of  the  custom  profile.  The 
results  are  shown  in  the  Chapter  VII  on  Data  Analysis. 

2.  Atmosphere  Effects  Module 

This  module  calculates  ray  refraction,  atmospheric  distortions  for  both  the  target 
and  the  scene,  transmissivity  and  path  radiance.  After  the  execution  of  the  ray  refraction 
calculation  (description  of  the  rays  that  arrive  at  the  sensor  or  target),  a  number  of  outputs 
are  available. 

The  Maximum  Intervision  Range  (MIVR)  (absolute  limit  of  detection  for  a  given 
target  height,  with  respect  only  to  ray  refraction  and  the  earth’s  curvature)  with  and 
without  refraction  (geometrical  horizon),  geometric  results  for  rays  to  ocean,  cloud  and 
target  are  some  of  the  outputs  of  that  module.  Elevation  angles  to  target,  sensor,  ocean 
intercept  and  cloud  intercept  are  calculated. 

Atmospheric  effects  such  as  inferior  and  superior  mirages  are  treated.  Effects  of 
atmospheric  turbulence  (blurring,  scintillation)  on  the  target  signature  based  on  the 
refractivity  structure  constant  are  also  calculated.  In  addition  atmosphere  target 
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magnification  and  demagnification  outputs  are  also  produced.  Moreover,  if  we  run  this 
module  in  the  "with  radiance"  mode,  additional  outputs  of  path  radiances  and 
transmissivities  from  sensor-to-target  and  sensor-to-background  are  produced.  For  those 
calculations  IRTOOL  uses  MODTRAN2  in  LOWTRAN  mode  with  a  selectable 
frequency  increment.  The  user  has  also  the  option  to  turn  the  refraction  effects  off.  In  this 
case  all  the  calculations  are  based  upon  the  straight  ray  geometry. 

Something  that  is  very  useful  to  the  user  is  the  ability  of  IRTOOL  to  manage  more 
than  one  passband.  For  each  of  them  the  module  calculates  the  values  for  each  target 
solution.  This  is  useful  for  the  prediction  of  other  band-dependent  quantities  such  as  path 
radiance  and  atmospheric  coherence  length. 

It  is  not  possible  to  run  the  atmosphere  effects  module  alone,  because  the  outputs 
of  the  Atmosphere  Profiler  Module  and  the  Engagement  Module  are  required.  After  the 
execution  of  those  two  modules  we  must  enter  a  specific  group  of  input  parameters. 
Those  parameters  are  shown  in  Table  12.  Parameters  which  are  shown  with  an  asterisk 
(*)  must  be  used  when  we  execute  the  Atmospheric  Effects  Module  in  the  "with 
radiance"  mode. 

3.  Engagement  Module 

The  Engagement  Module  produces  a  file  that  describes  the  motion  of  the  target 
relative  to  the  sensor  as  a  function  of  time.  The  user  can  specify  three  different  types  of 
trajectories.  In  the  first  mode  the  target  is  on  a  collision  course  with  the  sensor  and  moves 
with  constant  speed,  angle  of  attack  and  at  constant  altitude.  This  mode  has  two  user- 
selected  options.  In  the  first,  the  target  direction  of  motion  always  points  at  the  platform 
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position,  and  in  the  second  the  target  moves  to  the  predicted  position  of  target  at  the  time 
of  collision.  In  the  other  two  trajectory  modes,  the  user  must  specify  the  target  trajectory 
by  providing  it  in  an  appropriate  format.  In  the  first  one  we  have  to  specify  the  ground 
range  and  bearing  from  the  sensor  to  the  target  as  a  function  of  time,  and  in  the  second 
one,  the  user  must  specify  the  target  velocity  as  a  function  of  time  and  an  initial  position. 


PARAMETERS  GROUR  USED 

PARAMETERS 

Sensor  Type 

Number  of  Passbands  (min.  and  max. 
wavelength) 

Aperture  Diameter 

Sensor  Sampling* 

Sensor  Operation 

Sensor  Altitude 

Sensor  Look,  Azimuth 

Sensor  Look,  Elevation 

Environment 

Option  for  refraction  or  straight  ray  geometry 
Solar  Azimuth* 

Solar  Elevation* 

LOWTRAN  frequency  increment* 

Swell  rms  Amplitude 

Cloud  Base  Altitude* 

Cloud  Top  Altitude* 

OUTPUTS  NEEDED  FROM  OTHER 
MODULES 

Atmosphere  Profiler 

Engagement  (used  if  engagement  output  is 
available). 

Table  12.  Required  inputs  for  Atmosphere  El 

'fects  Module  in  IRTOOL. 

The  Engagement  Module  can  be  executed  alone  because  it  does  not  require  the 
outputs  from  other  modules.  The  input  parameters  required  for  the  successful  execution 
of  this  module  are  shown  in  the  Table  13.  The  trajectory  number  of  points  allows  us  to 
edit  the  number  of  points  calculated  along  the  target  trajectory.  This  number  affects  the 
execution  time  of  the  atmosphere  and  scintillation  modules.  Sometimes  the  user  under 
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certain  atmospheric  conditions  must  increase  this  number  in  order  to  obtain  reliable 
results.  From  [Ref.  2]  this  number  must  be  over  100  if  we  want  to  have  adequate  results. 


PARAMETERS  GROUR  USED 

PARAMETERS 

Sensor  Operation 

Sensor  Altitude 

Sensor  Look,  Azimuth 

Sensor  Look,  Elevation 

Sensor  Platform  Speed 

Sensor  Platform  Heading 

Target  Trajectory 

Trajectory  Mode 

Target  Maximum  Range 

Target  Minimum  Range 

Target  Altitude 

Target  Speed 

Target  Angle  of  Attack 

Initial  Bearing  to  Target 

Trajectory  #  of  points 

Table  13.  Required  inputs  for  the  Engagement  Module  in  IRTOOL. 


D.  RUNNING  OPTIONS 

The  Run  menu  allows  us  to  select  sequences  of  operations  (each  operation  is 
performed  by  an  entity  called  a  module),  referred  to  as  run  script.  Each  run  script 
corresponds  to  a  specific  set  of  input  parameters  in  memory  and  the  selected  modules.  A 
job  script  contains  one  or  more  run  scripts  that  run  sequentially  in  the  background. 

Under  the  run  menu  the  user  can  select  between  three  commands.  The  Default  Run 
Script,  the  Advanced  Run  Script  and  the  Any  Order  Run  Script. 

In  the  Default  Run  Script,  which  is  output  oriented,  the  user  selects  an  output  and 
IRTOOL  determines  what  calculations  need  to  be  performed.  A  list  of  modules  runs  with 
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the  present  set  of  input  parameters.  The  advanced  and  any  order  run  scripts  are  based  on 
modules,  where  each  module  is  responsible  for  doing  a  different  set  of  calculations.  The 
most  common  use  for  those  run  scripts  is  to  perform  data  analysis  or  processing  [Ref.  16]. 
In  the  advanced  run  script  the  user  chooses  what  modules  are  run.  The  modules  will  run 
in  the  order  that  they  appear.  In  the  any  order  run  script  which  also  contains  a  series  of 
modules  selected  by  the  user,  the  modules  can  be  run  in  any  order  determined  by  the  user. 
For  the  last  two  cases  we  must  have  in  mind  the  dependencies  between  the  different 
modules. 

For  this  research  we  used  mostly  the  default  run  script  “Engagement  (geometric), 
Atmosphere  and  Ray  Refraction”.  In  some  cases  we  used  the  Advanced  run  script, 
especially  when  we  had  to  use  a  LOWTRAN  file  created  by  the  user  (custom 
atmosphere).  Table  14  shows  which  default  run  script  we  should  use  in  order  to  get  the 
desired  output.  The  input  parameters  that  we  used  for  the  default  run  script  “Engagement 
(geometric),  Atmosphere  and  Ray  Refraction”  are  shown  in  the  Tables  11,12  and  13. 

E.  IRTOOL  CONVENIENCES 

IRTOOL  contains  user  friendly  interactive  IDL  programs  for  constructing  inputs 
needed  by  it  or  for  analysis  of  its  results.  These  conveniences  are  found  under  the 
Toolbox  menu.  The  IDL  routines  are  used  for  the  graphs  and  image  display  and  printing, 
for  editing  the  surface  layer  generated  by  the  Atmosphere  module,  for  converting 
measured  profiles  of  temperature,  pressure  and  relative  humidity  to  a  LOWTRAN  file 
usable  by  IRTOOL,  for  adding  our  own  interactive  applications  and  other  applications. 
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Desired  Output 

Input  Parameter  Groups 

Default  Run  Script 

Sensor  characterization 
(transfer  functions,  blip  noise, 
point  spread  function) 

Sensor 

Sensor 

Estimate  of  system 

Performance 

Sensor 

Sensor  operation 

Target 

Trajectory 

Environment 

Scene  Simulator 

System  Performance 

Algorithm  Control 

Estimated  Signal  to  Noise 

Ratio 

Marine  Boundary  Layer 
profile 

Environment 

Atmosphere  Profiler 

Maximum  Intervision  Range, 
boundary  profiles  and  horizon 
location 

Sensor 

Sensor  operation 

Environment 

Algorithm  Control 

Atmosphere  and  Ray 

Refraction  (horizon  location) 

Location  of  target  with 
refraction 

Sensor 

Sensor  operation 

Target 

Trajectory 

Environment 

Algorithm  Control 

Engagement  (geometric), 
Atmosphere  and  Ray 

Refraction 

Simulated  images  of 
background  without  sensor 
noise,  jitter,  or  blur 

Sensor 

Sensor  operation 

Environment 

Scene  Simulator 

Algorithm  Control 

Simulated  Scene-w/o  Sensor 
Effects 

Simulated  images  of 
background  with  sensor  noise, 
jitter,  and  blur 

Sensor 

Sensor  operation 

Environment 

Scene  Simulator 

Sensor  Effects 

Algorithm  Control 

Simulated  Scene-with  Sensor 
Effects 

Images  of  target  and  simulated 
background  without  sensor 
noise,  jitter  or  blur 

Sensor 

Sensor  operation 

Target 

Trajectory 

Environment 

Scene  Simulator 

Target  Injection 

Algorithm  Control 

Simulated  Scene-w/o  Sensor 
Effects,  with  target 
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Images  of  target  and  simulated 
background  with  sensor  noise, 
jitter  and  blur 

Sensor 

Sensor  operation 

Target 

Trajectory 

Environment 

Scene  Simulator 

Target  Injection 

Sensor  Effects 

Algorithm  Control 

Simulated  Scene-with  Sensor 
Effects,  with  target 

HDF  data  from  NATO  data 

IRAMMP  Data  Format 
Conversion 

Data  Analysis-IRAMMP  Data 
Format  Converter 

Reoriented,  subpatched  or 
processed  HDF  data 

Data  Preprocessor 

Data  Analysis-Preprocessor 

Mean,  variance,  correlation 
and  probability  densities  from 
HDF  data 

Data  Statistics 

•Data  Analysis-Statistics 

PSDs  and  cross  spectral 
densities  from  HDF  data 

Data  Spectral 

Data  Analysis-Spectral 

Table  14.  Inputs  and  Outputs  for  Default  Run  Scripts  in  IRTOOL  [  From  Ref.  16]. 


F.  OUTPUTS 

IRTOOL’ s  modules  produce  either  line  plot  or  image  data  files.  The  line  plot  files 
are  ASCII  files  with  the  extension  “.pit”.  They  are  readable  and  can  be  displayed.  The 
image  data  files  are  in  HDF  format  and  have  the  file  name  with  the  extension  “.img". 
Except  those  two  files,  IRTOOL  generates  error  (.err  extension)  and  log  (.log)  files.  The 
error  files  can  help  the  user  to  locate  a  problem  that  caused  the  module  to  terminate  its 
calculations.  The  log  files  contain  information  about  what  the  module  did. 

G.  NPS  PROFILER 

The  NPS  profiler  is  a  bulk  aerodynamic  Marine  Boundary  Layer  (MBL)  program 
which  uses  Monin-Obkhov  similarity  theory.  It  was  created  by  the  Meteorological 
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Department  of  the  NPS  and  was  provided  by  Prof.  K.Davidson.  The  program  is  written  in 
MATLAB  programming  language. 

The  meteorological  parameters  required  by  this  program  are  air  temperature,  sea 
temperature,  atmospheric  pressure,  wind  speed,  relative  humidity  and  the  heights  where 
these  measurements  were  made.  The  output  is  the  vertical  profiles  of  the  above 
parameters. 

We  used  the  NPS  atmosphere  profiler  because  this  code  requires  the  measurement 
height  for  each  input  parameter  in  contrast  with  IRTOOL,  that  requires  only  one  reference 
height,  assuming  that  all  measurements  were  taken  at  the  same  height.  But  as  we  will  see 
from  the  database  of  the  MAPTIP  trial  (Table  17)  for  some  datasets  this  was  not  the  case. 
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V.  MODEL  ALGORITHMS  FOR  REFRACTIVITY 

A.  INTRODUCTION 

Variations  of  pressure,  air  temperature  and  partial  pressure  of  the  various  gases  in 
the  atmosphere,  with  most  important  the  water  vapor  due  to  its  variation  with  height, 
cause  the  refraction  of  the  EM  waves.  In  order  to  model  the  refractivity  effects  in  the 
atmosphere  a  method  is  required  to  determine  the  vertical  profiles  of  temperature,  water 
vapor  pressure  (or  relative  humidity)  and  total  pressure. 


B.  IRBLEM 

The  L(W)WKD  is  a  marine  boundary  layer  (MBL)  model  which  calculates  the 
above  mentioned  profiles  based  on  the  similarity  (scaling)  theory  of  Monin  and  Obukhov 
[Ref.  17].  This  model  requires  measurements  of  the  water  temperature,  the  air 
temperature,  the  water  vapor  pressure  (or  relative  humidity),  the  total  pressure  and  wind 
speed.  The  fact  that  the  above  quantities  vary  principally  with  height,  allow  us  to  express 
the  refractive  index  as  a  function  only  of  height,  neglecting  any  horizontal  variations. 

The  refractivity  profile  and  its  first  derivative  can  be  determined  by  the  following 
equations  assuming  that  the  atmosphere  has  only  two  components,  a  dry  air  contribution 
and  a  water  vapor  contribution.  From  [Ref.  1 7,  Ref.  1 8] 

n(A,  /z)=l+R(T,  h)  xl(T6  (23) 
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The  first  derivatives  of  the  above  equations  are  given  by: 


ri(A,  h)= 


R’(A,  h)= 


^^-  =  R'(A,h)x  10"6 
dh 

dR(A, h)  (AK_ (h)  +  BV; (h)  -  R(A,  h)T\h)) 

dh  T{h) 


(25) 

(26) 


where  n  is  the  index  of  refraction 

ri  is  the  vertical  gradient  of  the  index  of  refraction  (dimensionless) 
R  is  the  refractivity  (dimensionless) 

R’  is  the  vertical  gradient  of  the  refractivity  in  units/meter 
A  is  the  wavelength  in  micrometers 

h  is  the  height  above  the  mean  water  level  (MWL)  in  meters 
Pa  =P-  Vp  is  the  dry  air  atmospheric  pressure  (mbars) 

P  is  the  atmospheric  pressure  (mbars) 

V p  is  the  water  vapor  pressure  (mbars) 


T  is  the  air  temperature  ( 0  K) 

A  ( A )  is  a  coefficient  given  by  Edlen  [Ref.  1 8]  for  the  dry  air  contribution 
B  ( A,  T)  is  the  coefficient  for  the  water  vapor  contribution 

Pa(h)  is  the  gradient  of  the  dry  air  pressure  in  mbars/meter 

T  (h)  is  the  gradient  of  the  temperature  and 

V'  ( h )  is  the  gradient  of  the  water  vapor  pressure. 
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From  [Ref.  18]  the  coefficients  A  and  B  can  be  assumed  constants  for  the  visible,  near 
infrared,  mid  infrared,  and  far  infrared.  Table  15  gives  the  values  of  those  coefficients  for 
these  bands. 


Wavelength  Band 

Range  (  /mi) 

A  (K/mbar) 

B  (K/mbar) 

Visible 

0.5-1 

78.5 

67.0 

Near  Infrared 

1-3 

77.7 

65.0 

Mid  Infrared 

3-5 

77.5 

65.0 

Far  Infrared 

8-12 

77.5 

40.0 

RADAR 

<20  GHz 

77.6 

373400/T 

Table  15.  Values  for  the  coefficients  A,  B  (given  by  Edlen)  that  are  used  in  IRBLEM  for 
the  calculation  of  refractivity  [From  Ref.  18]. 


With  the  use  of  the  Table  15,  Eq.  (26)  for  the  four  wavebands  can  be  written  as: 


Visible: 


dR(A,h)  _  78.5  dP  1S.5P  dT  11.5  dVp  U.5Vp  dT 
dh  T  dh  T 2  dh  T  ~dh  +  T2  ~dh 


(27) 


Near  Infrared: 


dR(A,h)  _  77.7  dP  11. IP  dT  12.2  dVp  \22Vp  dT 

dh  T  dh  T2  dh  T  dh  +  T2  dh 


(28) 


Mid  Infrared: 


dRjAJt)  _  77.5  dP  77.5?  dT_  12.5  dV p  12.5 dT 
dh  T  dh  T2  dh  T  dh  +  T2  dh 


(29) 


Far  Infrared: 


dR(X,h)  _  77.5  dP  77.5?  dT_  37.5  dVp  37.5Ep  dT 
dh  T  dh  T2  dh  T  dh  +  T2  ~dh 


(30) 
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C.  IRTOOL 

The  Atmosphere  Profiler  Module  using  the  bulk  method  and  the  Monin-Obukhov 
similarity  theory  produces  the  profiles  of  the  index  of  refraction  both  in  the  surface  layer 
and  also  in  the  region  above  it  [Ref.  11].  Assuming  again  no  horizontal  variations,  the 
index  of  refraction  is  given  by  n  =  1  +  Rx  10"6 ,  where  R  is  the  refractivity  and  from 
[Ref.  11] 

R  =RC  +  Ra  (31) 

where  Rc  is  the  continuum  value  and 

Ra  is  the  contribution  from  anomalous  dispersion  by  infrared  resonances  such 
as 

H  2  O,  CO 2 ,  O  3  etc  with  the  more  important  being  the  contribution  of  water. 
According  to  the  same  reference  “In  the  atmospheric  H  2  0  transmission  windows  the 
anomalous  dispersion  can  be  neglected.”  Under  this  assumption  Eq.  3 1  becomes 

R*  Rc  (32) 

The  continuum  value  of  refractivity  is  given  by  [Ref.  19]  as 

Rc=Rd+K  (33) 

where  Rd  is  the  contribution  from  dry  air  and 
Rw  is  the  contribution  of  water  vapor. 

From  [Ref.  1 1]  and  [Ref.  19]  the  value  of  Rd  and  Rw  is  given  by: 
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R,  --  °'™m  *.D  +  P-337  -  0.01 57TW)  -10"  P(h)]  (34) 

T(h) 

n  „00  29498.1  255.4 

with  i?0  =64.328  + - 7  H - 7  and  (35) 

0  146- r2  41 -r2 

Rw  =-1.765  xl0'l8(l-0.0109r2)Q(A)  (36) 

where  P  is  the  total  atmospheric  pressure  in  Torr 
T  is  the  temperature  ( 0  K) 

R  0  is  the  value  of  refractivity  at  a  temperature  of  1 5 0  C  and  pressure  of  760 
Torr  (1  atm) 

A  is  the  wavelength  in  micrometers  and 
Q  is  the  concentration  of  water  vapor  in  molecules/cm3 . 


For  atmosphere  modeling  IRTOOL  assumes  that  the  slight  wavelength  dependence  in 
Eqs.  35,  36  and  also  the  small  temperature  and  pressure  dependence  within  the  square 
bracket  in  Eq.  34  can  be  neglected.  Under  those  assumptions  the  Eqs.  34,  35,  36 
simplified  further  and  can  be  written  as: 


R  0-3789 mR 

d  T(h )  0 


Rw  =-1.765  xl0~'s  0(h) 


(38) 


R0  =  272.599  (IRTOOL  uses  the  value  272.7)  (39) 


Finally  the  refractivity  is  given  by: 
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R(h)  =0.3789 ^ - R0  - 1 .765  x  1 0~18 Q(h)  =1 03.326 - 1 .765  x  1  O’18 Q(h)  (40) 
V  '  T{h)  0  T(h) 


Using  the  Avogadro  constant  and  the  fact  that  the  coefficients  in  Eqs.  34  through  40  have 
the  units  necessary  to  make  the  equations  dimensionless  we  can  rewrite  Eq.  40  as: 

R(h)  =  n.5^-59Q  (41) 

where  P  is  the  total  pressure  in  mbars 
T  is  the  temperature  in  0  K  and 
Q  is  the  water  vapor  concentration  in  kg/m3 . 


The  ideal  gas  law  relates  the  water  partial  pressure  V  to  concentration  Q  by: 


v=0 — ^ — : 

p  FM(H20 ) 


(42) 


where  V  is  the  water  vapor  partial  pressure  in  Nt/m ' 


R  =8.314 


J 


(mole){K) 


is  the  universal  gas  constant 


kg 

FM(H20 )  =  0.018 — —  is  the  formula  mass  of  the  water  and 
mole 


T  is  the  temperature  in  degrees  Kelvin. 


With  the  use  of  Eq.  42  we  can  write  Eq.  41  as: 
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(43) 


The  vertical  profiles  of  the  parameters  in  Eq.  44  are  calculated  from  the  Atmosphere 
Profiler  module  of  IRTOOL. 

Here,  we  must  point  out  that  for  the  calculation  of  reffactivity  IRTOOL  does  not 
take  into  account  the  wavelength.  That  means  that  IRTOOL  does  not  model  dispersion 
(i.e.,  different  index  of  refraction  for  different  bands).  So,  there  is  no  difference  between 
the  ray  refraction  output  at  visible,  MWIR  or  LWIR.  In  contrast  in  IRBLEM  the 
wavelength  is  taken  into  account  through  the  coefficients  A  and  B.  The  use  of  the 
passband  in  IRTOOL  input  parameters,  as  we  showed  in  Chapter  IV,  is  only  for  band 
dependent  parameters  as  for  example  the  coherence  length  (r0). 

Under  the  assumption  that  the  earth  may  be  represented  by  a  sphere,  Snell’s  law 
which  describes  the  ray  refraction  in  the  atmosphere  can  be  written  [Ref.  20]  as: 

— («rcos/?)=0  (45) 

ds 

where  n  is  the  refractive  index  at  a  point 
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r  is  the  distance  of  that  point  from  the  earth’s  center  and 

/?  is  the  angle  between  a  particular  ray  and  the  local  horizon  plane  (Figure  17). 
Using  polar  coordinates,  the  trajectory  of  the  ray  can  be  described  with  the  use  of  the 
following  trigonometric  relation 

„  d\nr  ' 

tan/?  =  — — —  (46) 

uir 


With  the  use  of  Eqs.  45  and  46  we  have  that 


Figure  17.  Ray  Geometry  used  in  the  ray-tracing  program  of  IRTOOL  [After  Ref.  20] 


In  order  to  calculate  the  ray  trajectory  in  the  form  r(0,/?o)  where  J30  is  the  value  of  J3  for 

a  specific  ray  when  6  =  0 ,  we  use  the  Runga-Kutta  numerical  integration  method  [Ref. 
21]  in  Eqs.  46  and  47.  After  the  calculation  of  the  ray  trajectory  we  face  the  problem  of  if 
and  where  a  refracted  ray  actually  intercepts  the  ocean.  The  ocean  surface  supports  waves 
and  in  order  to  take  them  into  account  in  the  calculation,  IRTOOL  uses  the  significant 
wave  height  Hs  or  a  nominal  0.2  meters  whichever  is  greater  [Ref.  11].  The  significant 

wave  height  is  taken  to  be  equal  to  four  times  the  rms  surface  elevation  and  is  related  to 
wind  speed  by  an  approximate  relation  [Ref.  1 1,  Ref.  22] 

Hs=  0.02  lw2  (48) 

Where  u  is  the  wind  speed  in  m/sec  at  10  m  height.  If  in  addition  we  have  a  swell  then 
Eq.  48  is  modified  and  becomes 

0.021  [u2  +  (swell)2]  (49) 

This  is  the  cut-off  height  above  the  sea  surface  that  will  stop  all  the  rays. 
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VI.  EXPERIMENTAL  ARRANGEMENT  AND  PROCEDURE 


A.  INTRODUCTION 

In  1993  and  specifically  from  October  18  to  November  3  the  MAPTIP  (Marine 
Aerosol  Properties  and  Thermal  Imager  Performance)  trial  was  held  on  the  coast  of  The 
Netherlands.  It  was  sponsored  by  RSG  8  of  NATO’s  panel  4  and  the  principal  goals  of 
the  trial  were  to  study  the  performance  of  imaging  systems  within  the  marine  boundary 
layer  (MBL)  and  also  the  concentrations  and  properties  of  marine  aerosols  [Ref.  18]. 
After  the  trial  some  preliminary  results  were  presented  during  the  first  workshop  in  Oslo, 
Norway  (May  94)  by  most  of  the  participating  countries.  The  participants  also  agreed  to 
create  a  number  of  working  groups  so  as  to  jointly  study  and  combine  their  results.  Each 
group  had  a  specific  job  to  do  and  one  of  them  was  made  responsible  for  analyzing  and 
combining  the  results  of  the  various  refraction  measurements.  The  coordinator  of  this 
group  that  was  named  the  “Refractive  Effects  in  the  Visible  and  IR”  was  chosen  to  be  Dr. 
Forand  from  DREV  (Canada)  [Ref.  14], 

B.  EXPERIMENTAL  TECHNIQUES 

During  the  MAPTIP  trial  there  were  used  two  different  experimental  methods.  In 
the  first,  a  number  of  observation  devices,  namely  visible  and  IR  cameras,  were  placed  at 
known  heights  at  a  beach  station  and  a  spatial  array  of  radiation  sources  was  placed  at 
different  known  heights  above  the  MWL  on  a  platform  located  at  a  specific  distance  from 
the  beach  station.  This  arrangement  was  used  for  the  comparison  of  the  observed  angles 
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between  the  various  sources  with  both  the  geometric  angles  between  the  sources  and 
those  predicted  by  the  models  that  are  to  be  compared.  In  this  way,  keeping  the  range 
constant,  the  refractive  effects  of  varying  the  height  of  the  observers  or  the  sources 
(Height  Scan)  could  be  studied. 

In  the  second  method,  the  same  beach  station  was  used  while  the  radiation  sources 
were  placed  on  a  moving  platform  such  as  a  ship.  In  this  way,  the  elevations  of  the 
observers  and  the  sources  were  kept  constant,  while  the  range  between  the  station  and  the 
ship  was  continuously  varied  (Range  Scan).  This  method  was  used  for  the  recording  of 
two  parameters.  The  first  is  the  Maximum  Intervisibility  Range  (MIVR)  which  is  the 
range  at  which  a  source  at  a  particular  height  disappears  below  the  horizon,  and  the 
second  is  the  Minimum  Mirage  Range  (MMR)  which  is  the  minimum  range  at  which  a 
mirage  (if  there  is  one)  of  a  particular  source  begins  to  be  observed. 

C.  THE  EXPERIMENTAL  SETUP 

Figure  18  shows  the  geographical  location  of  the  beach  station  at  Katwijk  from 

which  the  observed  data  were  taken  (WPB),  the  location  of  the  Meetpost  Noordwijk 
(MPN)  platform  (WPN),  the  various  waypoints  (labeled  from  A  to  G)  and  the  tracks 
(marked  by  the  solid  lines)  followed  by  the  Hr.  Ms.  Tydeman.  Figure  19  shows  the  beach 
station  at  Katwijk  and  the  positions  of  the  different  cameras  that  were  used  during  the 
experiment.  The  cameras  were  placed  at  different  heights  above  mean  water  level  and 
their  types  and  characteristics  are  summarized  in  Table  16.  The  stationary  platform 
(Figure  20)  was  fixed  at  a  specific  distance  (10.4  km)  from  the  beach  station  and  eight 
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500W  halogen  lamps  were  mounted  to  the  platform  structure  at  heights  from  3.57  m  to 
19.84  m  above  the  mean  water  level.  The  lamps  (targets)  were  also  spaced  horizontally  so 
as  to  minimize  any  possible  overlap  of  the  recorded  images  and  secondary  images. 
Finally,  for  the  second  technique  the  Hr.  Ms.  Tydeman  was  used  as  a  moving  platform 
and  the  layout  of  the  sources  is  shown  in  Figure  21.  Six  500W  halogen  lamps  were 
mounted  at  various  heights  from  3.34  m  to  20.7  m  above  the  water  level  at  the  stem  of  the 
ship.  The  ship  was  asked  to  steer  in  a  straight  course  away  from  the  shoreline  while  a 
series  of  visible/IR  cameras  located  on  the  beach  station  were  used  to  record  the 
refraction  effects. 

D.  THE  MAPTIP  REFRACTION  DATASET 

Thirty  events  were  observed  during  the  MAPTIP  experiment.  Thirteen  of  those, 
labeled  from  A  to  M,  referred  to  Hr.  Ms.  Tydeman  (ship  tracking  events)  and  the 
remaining  seventeen  referred  to  the  stationary  MPN  platform  (MPN  observation  events). 
The  large  dataset  consists  of  video  recordings  obtained  with  the  cameras  we  mentioned 
above.  The  data  starting  with  the  letter  “T”  referred  to  the  ship  and  those  which  start  with 
the  letter  “P”  referred  to  the  platform. 

The  weather  conditions  were  also  obtained  using  a  weighted  average  from 
measurements  taken  by  different  participants  in  the  trial.  Due  to  the  large  dataset  it  was 
decided  that  only  a  portion  of  the  events  would  be  studied  in  great  detail.  This  subset 
consists  of  six  sets  of  data.  Each  set  of  data  consists  of  a  ship  tracking  event  and  two 
MPN  observation  events  taken  before  and  after  the  ship  event.  The  weather  conditions 
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and  especially  the  water  temperature,  air  temperature,  relative  humidity,  air  pressure, 
wind  speed  and  direction,  water  level  with  respect  to  MWL  (tide  height),  U/3  wave 
height,  solar  radiation,  rain  rate  and  the  heights  of  the  sensors  at  which  those  conditions 
were  measured  are  given  in  Table  17.  From  this  table,  we  observe  that  the  ASDT  (Air  to 
Surface  Difference  Temperature)  is  negative.  That  means  that  the  MBL  (Maritime 
Boundary  Layer)  is  unstable  and  sub-refractive  and  so  various  mirage  phenomena  are 
expected  to  occur.  Finally  Table  18  shows  the  cameras  that  were  used  for  each  of  the 
events  and  their  heights  above  the  mean  water  level. 


MAPTIP  Waypoints 


West  (km) 


Figure  18.  Map  of  the  geographic  location  of  the  MPN  and  the  various  waypoints  with 
respect  to  the  Beach  House  at  Katwijk  and  the  tracks  followed  by  the  Hr.  Ms.  Tydeman 
[From  Ref.  17]. 
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Figure  19.  Schematic  diagram  of  the  beach  station  showing  the  location  of  the  various 
cameras  during  the  MAPTIP  trial  [From  Ref.  17]. 
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Figure  20.  Schematic  diagram  of  the  MPN  platform  showing  the  placement  of  the  eight 
500  Watt  halogen  lamps  [From  Ref.  17]. 
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Figure  21 .  The  Hr.  Ms.  Tydeman  with  the  placement  of  the  six  halogen  lamps,  the 
exhaust  stack  and  the  stem  light  0.9  m  above  the  helicopter  deck  [From  Ref.  17]. 


Camera 

Country' 

Type 

Wavelength  (pm) 

Resolution 

(prad/pixel) 

VI 

Canada 

Sony  CCD 

>0.85 

6.23  ±2% 

V2 

Canada 

Sony  CCD 

>0.715 

9. 15  ±2% 

§■1 

Canada 

Mitsubishi 

3-5 

9.91  ±2% 

V3;G 

Germany  . 

CCD 

Visible 

11. 94  ±2% 

V4;F 

France 

Sony  CCD 

Visible 

161  ±5% 

IR2;  F 

France 

Castor 

8-12 

139±5% 

IR3;  F 

France 

Mitsubishi 

3-5 

468  ±5% 

Table  16.  1 

'ypes  and  characteristics  of  the  visible  and  infrared  cameras  that  were  used 

during  the  MAPTIP  trial  [After  Ref.  17]. 
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SET  1 

SET  2 

SET  3 

Parameters 

P-F2 

T-B* 

P-F3 

P-13 

T-D* 

P-14 

P-J2 

T-E 

P-J5 

Date  (dd-mm-yy) 

20-10-93 

20-10-93 

20-10-93 

22-10-93 

22-10-93 

22-10-93 

22-10-93 

22-10-93 

22-10-93 

Start  (hh:mm) 

7:29 

7:25 

8:12 

14:46 

14:45 

15:33 

20:22 

20:45 

21:59 

Stop  (hh:mm) 

7:34 

8:25 

8:19 

14:51 

15:40 

15:36 

20:24 

21:55 

22:01 

Water  Temp.  (C) 

12.7 

12.5 

13.1 

12 

12.4 

12.3 

12 

12.4 

12.1 

!■■■ 

Ht.  of  Sensor  1  (m) 

12 

12 

12 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

Air  Temp.  (C) 

8.3 

8.4 

8.6 

9.2 

9.1 

8.9 

7.4 

7.4 

7 

Rel.  Humidity  (%) 

65 

65 

64 

69 

72 

73 

79 

76 

73 

Air  Pressure  (mbar) 

1026 

1026 

1026 

1026" 

1025 

1026 

1028 

1027 

1027 

Ht.  of  Sensor  2  (m) 

22.6 

22.6 

22.6 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

Wind  Speed  (m/s) 

4.5 

5 

6.2 

10" 

9.8 

10 

9 

9.1 

9.5 

Wind  Direction  (deg) 

287 

280 

210 

60 

50 

45 

62 

60 

60 

Water  Level  (m) 

0.44 

0.33 

0.13 

-1.02 

-1.02 

-1.02 

0.5 

0.32 

0.13 

Wave  Height  (m) 

0.18 

0.17 

0.16 

0.98 

0.98 

1 

0.94 

0.91 

0.86 

Solar  Radiat.  (W/m2) 

80 

70 

70 

42 

30 

15 

0 

0 

0 

Rain  Rate  (mm/lOmin) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ASTD(C) 

-4.4 

-4.1 

-4.5 

-2.8 

-3.3 

-3.4 

-4.6 

-5 

-5.1 

SET  4 

SET  5 

SET  6  | 

Parameters 

P-Ml 

T-F 

P-M2 

P-M9 

T-G 

P-M10 

P-P2 

T-K 

P-P3 

Date  (dd-mm-yy) 

25-10-93 

25-10-93 

25-10-93 

25-10-93 

25-10-93 

25-10-93 

28-10-93 

28-10-93 

28-10-93 

Start  (hh:mm) 

3:17 

3:35 

4:50 

14:53" 

15:40 

16:34 

8:08 

8:35 

9:18 

Stop  (hh:mm) 

3:18 

4:40 

4:53 

14:55 

16:25 

16:37 

8:10 

9:35 

9:24 

Water  Temp.  (C) 

11.6 

11.6 

11.6 

11.6 

11.5 

11.5 

11.6 

11.7 

11.7 

naiiH 

■■■■ 

Ht.  of  Sensor  I(m) 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

Air  Temp.  (C) 

10.8 

11 

11 

10.2 

10.3 

10 

8.9 

9.2 

9 

Rel.  Humidity  (%) 

72 

73 

74 

76 

77 

77 

80 

80 

78 

Air  Pressure  (mbar) 

1028 

1030 

1029.5 

1031 

1031 

1031 

1027 

1027 

1027 

Ht.  of  Sensor  2  (m) 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

Wind  Speed  (m/s) 

8 

7.3 

7.5 

7 

8.2 

7 

4.8 

4.1 

4.1 

Wind  Direction  (deg) 

30 

25 

30 

40 

30 

35 

70 

82 

80 

Water  Level  (m) 

-0.47 

-0.51 

-0.67 

-0.78 

-0.78 

-0.75 

-0.79 

-0.79 

Wave  Height  (m) 

0.9 

0.92 

0.97 

0.91 

0.93 

0.92 

0.51 

0.51 

0.51 

Solar  Radiat.  (W/m2) 

0 

0 

0 

20 

5 

0 

50 

75 

90 

Rain  Rate  (mm/lOmin) 

0 

0 

0 

0 

0 

0 

0 

0 

ASTD  (C) 

-0.8 

-0.6 

-0.6 

-1.4 

-1.2 

-1.5 

-2.7 

-2.7 

tmam 

■■■ 

■■■ 

asm 

■■■1 

^■n 

■■■ 

mam 

mam\ 

Note:  *  -WPC  course 

1 

i 

[ 

1 

Table  17.  The  dataset  (events  and  their  associated  meteorological  conditions)  that  was 
used  for  the  purposes  of  analysis  [From  Ref.  17]. 
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Events 


Visible 


1  Visible  2 


IR1  (3-5) 


Camera 


P2010F2 

P2010F3 

T2010B 

P2210I3 

P2210I4 

T2210D 

P2210J2 

P2210J5 

T2210E 

P2510M1 

P2510M2 

T2510F 

P2510M9 

P2510M10 

T2510G 

P2810P2 

P2810P3 

T2810K 


(Canada) 
7.04  m 
7.04  m 


7.04  m 
7.04  m 


7.04  m 
7.04  m 
7.04  m 
7.04  m 
7.04  m 
7.04  M 
7.04  m 
7.04  m 
7.04  m 
7.04  m 


7.04  m 


Camera 

(Canada) 


6.44  m 


6.39  m 
8.91  m 
8.91  m 


3.85  m 
3.85  m 
8.83  m 
3.85  m 
3.85  m 
3.85  m 
3.85  m 
3.85  m 
3.85  m 


Camera 

(Canada) 


6.89  m 


7.04  m 


7.04  m 
7.04  m 
7.04  m 
7.04  m 
7.04  m 
7.04  m 


7.04  m 
7.04  m 


Visible  3 
Camera 
(Germany) 


Visible  4 
Camera 
(France) 


IR2  (8-12) 

Camera 

(France) 


20  m 


20  m 


20  m 


20  m 


6.63  m 


6.63  m 


6.63  m 


6.79  m 


6.79  m 


6.79  m 


IR3  (3-5) 

Camera 

(France) 


6.71  m 


6.71  m 


6.71  m 


Table  1 8.  The  cameras  that  were  used  for  each  of  the  events  during  the  MAPTIP  trial  and 
their  heights  above  the  mean  water  level  [After  Ref.  18]. 


VII. 


DATA  ANALYSIS 


A.  INTRODUCTION 

A  comparison  of  IRBLEM  with  the  French  model  PIRAM,  using  the  same 
database  shown  in  Table  17,  was  described  in  1997.  The  results  were  very  satisfactory  for 
both  but  better  for  IRBLEM  [Ref.  17]. 

In  1998,  an  effort  was  made  to  compare  IROOL,  IRBLEM  and  PIRAM  under  the 
same  set  of  input  conditions  [Ref.  2],  The  ship  tracking  events  of  the  MAPTIP 
experiment  were  used  and  the  analysis  showed  that  IRTOOL  produces  MMRs  at  greater 
ranges  than  both  the  observed  ones  and  those  predicted  from  the  other  two  models.  As  is 
known,  for  a  target  at  a  constant  altitude  under  subrefractive  conditions  a  mirage  appears 
only  when  this  target  is  in  a  specific  range  zone  limited  by  the  MMR  and  MIVR.  When 
the  MMR  moves  towards  greater  ranges  the  range  zone  gets  narrower.  As  a  result,  some 
of  the  targets  that  were  in  the  range  zone  before,  may  be  now  out  of  it  and  mirages  will 
not  be  predicted  for  them.  So,  we  expect  that  IRTOOL  may  not  predict  mirage  for  some 
targets  for  which  mirages  were  found  experimentally.  As  mentioned  in  Chapter  I,  Forand 
at  DREV  applied  IRTOOL  to  data  obtained  off  the  Netherlands  coast  in  the  MAPTIP 
measurements  in  1993  and  concluded  that  IRTOOL  did  not  predict  mirage  in  cases  in 
which  they  were  found  experimentally  (L.  Forand,  private  communication).  According  to 
the  discussion  above,  this  opinion  is  in  agreement  with  the  results  that  were  obtained 
during  the  analysis  in  [Ref.  2], 
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B.  DATA  ANALYSIS 

In  this  chapter  we  will  compare  the  IRTOOL  and  IRBLEM  models  and  also  we 
will  examine  again  the  capability  of  IRTOOL  for  producing  mirage  effects,  using  the 
database  shown  in  Table  17. 

The  appearance  of  the  mirage  effects  depends  on  refractivity.  For  that  reason  we 
will  examine  the  profile  of  refractivity  as  it  is  produced  from  the  two  models.  Also  we 
will  examine  the  profiles  of  temperature,  humidity  and  pressure  because  refractivity 
depends  on  those  profiles.  We  will  discuss  any  differences  in  the  results  and  we  will 
make  suggestions  for  IRTOOL  such  as  to  improve  its  results  and  to  get  the  best  fit 
agreement  with  the  experimental  data. 

In  order  to  be  consistent  with  the  previous  two  analyses  we  will  use  the  IRBLEM 
program  as  it  was  used  in  [Ref.  18]  and  instead  of  extrapolating  the  results  from  that 
reference  as  was  done  in  [Ref.  2]  we  will  run  the  IRBLEM  code.  According  to  [Ref.  1 8] 
IRBLEM  and  PIRAM  were  employed  using  a  wave  height  equal  to  1.3  times  the 
measured  Hy  wave  height  (the  “significant  wave  height”)  and  also  equal  to  1.0  times  the 
measured  H  y.  Moreover,  in  the  analysis  both  modes  of  the  refraction  module  of 
IRBLEM  were  used,  the  “wavy”  (L(W)WKD)  and  the  “non-wavy”  mode  (LWKD).  The 
conclusion  was  that  the  combination  of  a  wave  height  equal  to  1 .3  times  the  measured 
wave  height  with  the  non-wavy  mode  produces  results  closer  to  the  observed  data  [Ref. 
17].  So,  in  our  analysis  we  will  use  those  conditions.  After  the  completion  of  the  data 
analysis  it  will  possible  to  compare  all  the  models  (IRTOOL,  IRBLEM,  and  PIRAM) 
using  the  results  from  this  thesis  and  the  results  from  [Ref.  18].  The  results  of  our  analysis 
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are  presented  here  referenced  by  the  “Event  index”,  referring  to  the  measurements  of  the 
MAPTIP  experiment. 

Before  we  proceed,  it  is  very  important  to  note  that  the  wave  height  in  Table  17  is 
the  measured  H/3  wave  height  and  not  the  rms  wave  height  [Ref.  18].  For  that  analysis 
all  the  heights  we  will  refer  to  are  above  the  water  level  (i.e.,  the  height  of  the  tide  is 
taken  into  account).  The  height  of  the  tide  is  also  shown  in  the  Table  17  (water  level). 
Events  starting  with  the  letter  P  refer  to  the  stationary  platform  and  those  with  the  letter  T 
refer  to  the  moving  platform  (ship). 

1.  P2010F2 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (VI)  camera  positioned  at  height  h  v  =  6.6  m  above  the  water 

level.  The  comparison  of  the  computations  for  this  event  is  shown  in  Figure  22,  where  the 
solid  line  with  the  circles  represents  a  non-refractive  atmosphere,  the  squares  represents 
the  observed  data,  the  rhombus  represents  the  prediction  from  the  IRBLEM  code  and  the 
triangles  the  data  predicted  from  the  IRTOOL  code.  As  we  can  see,  IRBLEM  predicts 
mirages  for  the  same  targets  as  the  experimental  data,  while  IRTOOL  does  not  predict  the 
mirage  for  the  target  at  7.19  m  and  also  does  not  detect  the  target  at  3.13  m.  We  can  also 
observe  that  any  target  below  «  3  m  will  not  be  detected  from  IRBLEM.  This  height 
increases  to  «3.5  m  for  IRTOOL.  For  a  non-refractive  atmosphere  any  target  below  «0.5 
m  should  be  below  the  horizon.  Finally,  the  shape  of  the  curves  that  both  models  produce 
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deviates  from  the  observations  in  the  region  of  their  cusps.  It  is  also  clear  from  Figure  22 
that  the  results  predicted  from  IRBLEM  are  much  better  than  from  IRTOOL. 

As  we  have  already  mentioned  in  Chapter  IV,  IRTOOL  does  not  use  the 
measurement  of  pressure  as  an  input  parameter  even  though  the  pressure  is  used  in  the 
calculation  of  refractivity  (Eqs.  43,  44).  The  same  model  also  assumes  that  all  the 
measurements  of  temperature,  humidity,  pressure  and  wind  speed  are  made  at  the  same 
height.  From  Table  17  the  air  temperature  and  humidity  were  measured  at  12  m  while  the 
anemometer  was  placed  at  a  height  of  22.6  m.  For  those  two  reasons  we  used  the 
atmosphere  profiler  of  IRTOOL  to  create  a  custom  profile.  This  profile  was  created  using 
the  NPS  profiler  in  IRTOOL,  in  which  the  pressure  is  an  input  and  also  requires  the 
measurement  height  for  each  parameter.  In  Figures  23A,  23B,  23C  we  see  the  profiles  of 
humidity,  temperature  and  pressure  generated  from  the  atmosphere  profiler  module  of 
IRTOOL  using  the  midlatitude  winter  atmosphere  (symbol  “x”)  and  the  custom 
atmosphere  (symbol  “o”).  The  generation  of  this  custom  atmosphere  was  made  as 
follows: 

First,  we  run  the  NPS  profiler  with  the  same  input  parameters  to  create  the  profiles  of 
temperature,  relative  humidity,  pressure  and  wind  speed.  These  profiles  are  then  used  in 
conjunction  with  the  Atmosphere  Shaper  under  the  Toolbox  menu  of  IRTOOL  to  produce 
the  custom  atmosphere.  In  Figures  23A,  23B,  23C  are  also  shown  the  profiles  as  they 
were  calculated  from  the  LWKD  refraction  model  of  IRBLEM  (solid  line).  In  Figure  23  A 
(humidity  profile)  we  see  no  difference  between  the  profiles  except  that  at  the  surface  the 
relative  humidity  calculated  from  IRBLEM  is  108%  which  is  a  non-physical  result.  On 
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the  other  hand  in  IRTOOL  the  relative  humidity  cannot  exceed  100%.  In  Figure  23B 
(temperature  profile)  the  profiles  of  the  midlatitude  winter  atmosphere  and  custom 
atmosphere  are  identical  but  there  is  a  slight  difference  from  the  profile  predicted  from 
IRBLEM,  especially  in  the  first  few  meters  above  the  water  surface.  Here  we  must  note 
that  IRTOOL  calculates  the  right  temperature  of  the  sea  surface  (12.7°C)  in  contrast  with 
IRBLEM,  which  gives  a  temperature  of  11.58°C.  In  Figure  23C  (pressure  profile)  we 
observe  that  the  profile  of  pressure  predicted  from  IRBLEM  and  that  which  is  predicted 
from  IRTOOL  using  the  custom  atmosphere  are  exactly  the  same,  unlike  the  IRTOOL 
with  the  midlatitude  winter  atmosphere,  which  predicts  different  values  for  the  pressure. 
There  is  a  difference  of  9.5  mbar  or  7.12  Torrs.  In  spite  of  this  both  models  give  the  same 
pressure  gradient.  That  means  the  error  in  using  the  midlatitude  winter  atmosphere 
instead  of  the  custom  atmosphere  for  the  calculation  of  refractivity  will  be  less  than 
0.93%  in  accordance  with  Eq.  43  and  Figure  23C.  From  Figure  23D  we  observe  also  a 
small  difference  in  the  wind  speed  profiles  between  our  models.  Here  we  must  notice  that 
the  wind  speed  at  the  water  surface  calculated  from  IRBLEM  is  0.72  m/sec.  On  the  other 
hand  IRTOOL  always  assumes  that  the  wind  speed  at  the  sea  surface  is  0  m  /sec.  This  is 
another  difference  between  the  two  models. 

As  we  have  already  mentioned,  the  refraction  (ray  bending)  depends  on  the 
gradient  of  refractivity  and  not  on  the  absolute  magnitude  of  the  refractivity  itself.  That 
means  that  we  are  interested  in  the  final  result,  which  is  the  profile  of  refractivity  or  its 
gradient.  In  Figure  23E  is  shown  the  profile  of  refractivity.  We  see  that  for  the  same 
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height,  the  custom  atmosphere  and  the  midlatitude  winter  atmosphere  lead  to  different 
values  for  the  refractivity  but  from  Figure  23F  the  gradient  of  reffactivity  for  those  two 
atmospheres  is  the  same.  That  means  the  error  in  pressure  is  neglible  in  the  calculation  of 
reffactivity  and  also  the  restriction  to  use  only  one  reference  height  does  not  affect  at  all 
the  final  result  which  is  the  gradient  of  reffactivity,  at  least  for  that  event.  From  the  same 
Figure  23F,  we  see  a  slight  difference  in  the  first  few  meters  between  IRTOOL  and 
IRBLEM.  We  observe  that  the  first  derivative  of  refractivity  for  IRTOOL  is  greater  than 
that  of  IRBLEM.  That  means,  the  path  curvature  of  the  rays  in  IRTOOL  for  the  first  few 
meters  is  stronger  or  the  radius  of  curvature  of  the  ray  is  smaller  (Eqs.  6,  7).  So,  the 
refraction  effects  in  IRTOOL  are  stronger  for  the  first  few  meters.  So,  we  expect  that  the 
MMR  (Minimum  Mirage  Range)  of  IRTOOL  moves  toward  shorter  ranges.  Also,  this 
strength  of  reffaction  is  the  reason  that  IRTOOL  does  not  predict  an  image  of  the  target  at 
3.13  m.  This  can  be  understood  by  looking  in  Figure  9.  But  if  the  MMR  in  IRTOOL 
moves  to  shorter  ranges  than  IRBLEM  then  we  should  expect  a  mirage  at  least  for  the 
target  at  7.19  m.  The  reason  that  this  does  not  happen  must  lie  in  the  cut-off  height.  We 
recall  from  Chapter  V  that  another  significant  factor  that  affects  the  propagation  of  the 
low  level  rays  is  the  significant  wave  height  or  cut-off  height.  An  increased  wave  height 
would  move  the  predicted  MMR  to  longer  ranges.  The  significant  wave  height  that  was 
calculated  and  used  for  the  ray  tracing  program  in  IRTOOL  was  H5  =0.5 1  m.  On  the  other 

hand  in  the  ray  tracing  program  of  IRBLEM  the  height  of  H5  =0.1 7  m  was  used  according 
to  the  Eq.  19  with  H„,=1.3xHX  and  minabslayer=0.03  m.  The  difference  between  the 
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cut-off  heights  that  were  used  in  the  ray  tracing  programs  is  responsible  for  the  failure  of 
IRTOOL  to  predict  the  mirage  for  the  target  at  7.19  m.  In  order  to  make  sure  that  this  is 
true,  we  edit  the  FWH  of  IRBLEM  so  as  to  produce  the  cut-off  height  of  IRTOOL.  For 
that  reason  we  used  the  value  of  FWH=2.05.  The  reason  that  we  edit  the  cut-off  height  of 
IRBLEM  instead  of  IRTOOL  is  that  in  IRBLEM  if  we  change  the  FWH  we  do  not  change 
the  calculated  profiles.  On  the  other  hand  the  only  way  to  change  the  cut-off  height  in 
IRTOOL  is  to  use  a  wind  speed  different  than  the  measured  one.  But  doing  this  will 
change  the  profiles  and  so  we  will  not  be  able  to  compare  those  models  under  the  same 
set  of  input  conditions.  Previously,  we  said  that  the  restriction  of  IRTOOL  to  use  only  one 
reference  height  does  not  affect  at  all  the  calculation  of  the  reffactivity  gradient.  This  was 
true  at  least  for  that  event.  On  the  other  hand,  the  use  of  only  one  reference  height  will 
affect  the  calculation  of  the  cut-off  height  which  is  a  crucial  factor  for  the  ray-tracing. 
That  happens  because  the  wind  speed  at  10  m  height  is  used  for  the  calculation  of  the  cut¬ 
off  height  in  IRTOOL.  In  this  event  we  used  as  a  reference  height  the  height  of  12  m 
because  most  of  the  measurements  were  made  at  this  height.  But  the  anemometer  was  at 
the  height  of  22.6  m.  That  means  the  cut-off  height  that  was  calculated  was  not  correct.  If 
we  had  used  as  reference  height  the  height  of  22.6  m  then  the  extrapolated  wind  at  10  m 
would  be  lower,  and  as  a  result  the  cut-off  height  would  be  smaller  because  in  the 
calculation  of  the  significant  height  the  wind  speed  was  used  in  quadrature.  The  results 
with  the  edit  cut-off  height  of  IRBLEM  are  shown  in  Figure  24.  Now,  IRBLEM  does  not 
predict  the  mirage  for  the  target  at  7.19  m,  but  the  target  at  3.13  m  is  still  detectable. 
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Figure  25  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  06:45  to  08:25  GMT.  The  event  that  we  analyzed  took  place  from  07:29  to  07:34 
GMT  and  the  ASTD  that  was  used  was  ASTD  =  -4.4 0  C.  According  to  Figure  25  there  is 
a  variation  in  ASTD  of  approximately  2 0  C  from  10-20  km.  That  means  the  results  can  be 
further  improved  if  we  take  into  consideration  this  figure.  Using  a  lower  magnitude  of 
ASTD  the  strength  of  the  refraction  for  both  models  will  be  reduced,  and  the  results  will 

be  moved  closer  to  the  observed  data.  Figure  26  shows  the  results  for  an  ASTD  =  -3.9 0  C, 
and  using  the  same  cut-off  height  for  both  models.  The  target  at  3.13  m  is  not  detected 
from  IRTOOL,  but  this  time  IRBLEM  does  not  predict  the  mirage  of  the  target  at  4.58  m. 
Decreasing  the  refraction  even  more  by  using  an  ASTD  =  -3.7 0  C  (Figure  27)  we  observe 
that  IRTOOL  now  detects  the  target  at  3.13  m.  The  results  are  moving  even  closer  to  the 
observations  but  for  the  target  at  4.58  m  IRBLEM  does  not  predict  the  mirage.  With  an 
ASTD  =  -3.4°  C  (Figure  28)  there  is  very  good  agreement  between  the  predicted  and  the 
observed  data,  and  both  models  give  approximately  the  same  solutions.  It  is  important  to 
note  that  the  change  in  ASTD  affects  IRTOOL  more  than  IRBLEM  as  we  can  see  from 
Figures  24,  26,  27  and  28. 

2.  P2010F3 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (VI)  positioned  at  height  h5  =  6.91  m  above  the  water  level. 

From  Table  17  we  see  that  the  air  temperature  and  humidity  were  measured  at  12 
m  while  the  anemometer  was  placed  at  a  height  of  22.6  m.  For  the  reference  height  in 
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IRTOOL  we  used  12  m.  This  will  produce  an  error  less  than  3%  in  the  value  of  the  wind 
speed  at  the  height  of  the  22.6  m,  and  so  we  do  not  expect  the  results  to  change 
significantly.  But  in  order  to  minimize  any  error  we  generated  the  wind  speed  profile 
from  the  NPS  profiler.  The  wind  speed  that  was  calculated  at  12  m  was  6.03  m/sec,  and 
that  was  the  value  of  the  wind  speed  that  we  used  in  the  input  parameters  of  IRTOOL. 
The  value  of  the  wind  speed  at  12  m  as  it  was  predicted  from  the  NPS  profiler  is  in 
agreement  with  that  of  IRBLEM. 

As  we  can  see  in  Figure  29,  IRBLEM  gives  mirages  for  the  same  targets  as  the 
experimental  data,  when  IRTOOL  does  not  predict  the  mirages  for  the  targets  at  7.5  m 
and  4.89  m.  Also,  it  does  not  detect  the  target  at  3.44  m.  We  can  also  observe  that  any 
target  below  «  3  m  will  not  be  detected  from  IRBLEM.  This  height  increases  to  «3.5  m 
for  IRTOOL:  For  a  non-refractive  atmosphere  any  target  below  «1  m  should  be  below 
the  horizon.  Finally,  the  shape  of  the  curves  that  both  models  produce  does  not  follow  the 
data  in  the  region  of  their  cusps.  From  the  same  figure  it  is  also  clear  that  the  results 
predicted  from  IRBLEM  are  much  better  than  from  IRTOOL. 

In  Figures  30A,  30B,  30C,  30D,  30E,  30F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  as  they  were 
generated  from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (symbol  “x”)  and  the  custom  atmosphere  (symbol  “o”)  and  also  those  that 
were  generated  from  the  refraction  module  of  IRBLEM  (solid  line).  From  the  Figures 
30A  and  30B  we  observe  that  at  0  m  the  value  of  the  humidity  is  again  greater  than  100% 
and  the  temperature  1  °C  less  than  the  measured  one.  In  IRTOOL  for  the  same  height  the 
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relative  humidity  is  100%  and  the  water  temperature  13.1  °C  which  is  in  agreement  with 
the  value  given  in  Table  17.  From  Figure  30D  we  observe  a  small  difference  in  the  wind 
speed  profiles  and  also  that  the  wind  speed  at  the  water  surface  calculated  from  IRBLEM 
is  0.53  m/sec.  From  Figure  30F  we  observe  that  again  the  two  atmospheres  that  we  used 
in  IRTOOL  give  exactly  the  same  results  for  the  gradient  of  refractivity  (ray  bending). 
The  error  in  pressure  using  the  midlatitude  winter  atmosphere  (Figure  30C)  is  about  0.9% 
and  so  it  is  not  going  to  affect  the  profile  of  the  gradient  of  refractivity.  Also,  we  observe 
a  slight  difference  in  the  first  few  meters  between  IRTOOL  and  IRBLEM.  The  first 
derivative  of  refractivity  in  IRTOOL  is  greater  than  that  in  IRBLEM.  This  means  that  the 
path  curvature  of  the  rays  in  IRTOOL  for  the  first  few  meters  is  stronger  or  in  another 
words  the  radius  of  curvature  of  the  ray  is  smaller.  So,  we  expect  that  the  MMR 
(Minimum  Mirage  Range)  of  IRTOOL  moves  toward  shorter  ranges.  This  strength  of 
refraction  is  the  reason  that  IRTOOL  does  not  predict  the  target  at  3.44  m.  But  if  the 
MMR  in  IRTOOL  moves  to  shorter  ranges  than  IRBLEM  then  we  should  expect  a  mirage 
at  least  for  the  targets  at  4.89  m  and  7.5  m.  So,  we  edited  the  cut-off  height  in  IRBLEM 
(FWH=4.18)  for  the  same  reasons  that  we  mentioned  in  the  event  P2010F2.  The 
significant  wave  height  that  was  calculated  and  used  from  IRTOOL  was  Hs  =0.90  m  and 

for  IRBLEM  was  H  v=0.155  m  (actually  the  same  as  the  measured  wave  height).  Even 

though  we  used  a  lower  value  for  the  wind  speed  in  the  input  parameters  of  IRTOOL,  the 
calculated  significant  wave  height  is  still  much  greater  than  the  measured  one  (error 
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82%).  The  results  with  the  edit  cut-off  height  of  IRJBLEM  are  shown  in  Figure  3 1 .  Here 
only  IRBLEM  predicts  an  image  for  the  target  at  3.44  m,  which  is  still  detectable. 

Figure  25  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  06:45  to  08:25  GMT.  The  event  that  we  analyze  took  place  from  08:12  to  08:19 
GMT  and  the  ASTD  that  was  used  was  ASTD  =  -4.5  °C.  According  to  Figure  25  there  is 
a  variation  in  ASTD  approximately  2  0  C  from  10-20  km.  That  means  the  results  can  be 
further  improved  if  we  take  into  consideration  this  figure.  Using  a  lower  magnitude  of 
ASTD  the  strength  of  the  refraction  for  both  models  will  be  reduced,  moving  the  results 
closer  to  the  observed  data.  In  this  case  we  expect  that  for  a  specific  ASTD  the  IRTOOL 
will  detect  the  target  at  3.44  m,  but  we  are  not  sure  if  it  will  continue  to  predict  the  mirage 
of  the  target  at  4.19  m  (unless  if  we  run  the  IRTOOL  code  for  each  ASTD).  On  the  other 
hand  the  sure  thing  is  that  IRBLEM  will  continue  not  to  predict  the  mirage  of  any  target, 
except  maybe  the  target  at  3.44  m. 

3.  P2810P2 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (V2)  positioned  at  height  h5  =  4.6  m  above  the  water  level. 

From  Table  17  we  see  that  all  the  meteorological  measurements  were  made  at  the 
same  height  of  3.4  m.  That  means  we  do  not  have  to  use  the  NPS  profiler  to  get  the  value 
of  the  wind  speed  at  3.4  m.  As  we  can  see  in  Figure  32,  IRBLEM  gives  mirages  for  the 
same  targets  as  the  experimental  data  while  IRTOOL  does  not  predict  the  mirage  for  the 
target  at  5.77  m.  We  can  also  observe  that  any  target  below  «  3.8  m  will  not  be  detected 
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by  IRBLEM.  This  height  increases  to  «4  m  for  IRTOOL.  For  a  non-refractive 
atmosphere  any  target  below  « 1  m  should  be  below  the  horizon.  Finally,  the  shape  of  the 
curves  that  both  models  produce  does  not  follow  the  data  in  the  region  of  their  cusps. 
From  the  same  figure  it  is  also  clear,  that  the  IRBLEM  prediction  results  are  much  better 
than  IRTOOL. 

In  Figures  33A,  33B,  33C,  33D,  33E,  33F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  reffactivity  and  gradient  of  refractivity  as  they  were 
generated  from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (symbol  “x”)  and  the  custom  atmosphere  (symbol  “o”)  and  also  those  that 
were  generated  from  the  refraction  module  of  IRBLEM  (solid  line).  From  the  Figures 
33  A  and  33B  we  observe  that  at  0  m  the  value  of  the  humidity  is  again  greater  than  100% 
and  the  temperature  0.7  °C  less  than  the  measured  one.  In  IRTOOL  for  the  same  height 
the  relative  humidity  is  100%  and  the  water  temperature  11.6°C  which  is  in  agreement 
with  the  value  given  in  Table  17.  From  Figure  33  D  we  observe  a  small  difference  in  the 
wind  speed  profiles  and  also  that  the  wind  speed  at  the  sea  surface  calculated  from 
IRBLEM  is  0.7  m/sec.  From  Figure  33C  we  observe  that  the  error  in  the  second  term  of 
Eq.  44  using  the  midlatitude  winter  atmosphere  instead  the  custom  atmosphere  for  the 
calculation  of  the  pressure  profile  will  not  be  greater  than  0.90%.  This  is  shown  in  Figure 
33F.  We  observe  that  again  the  two  atmospheres  that  were  used  in  IRTOOL  give  the 
same  results  for  the  gradient  of  refractivity  (ray  bending).  Also,  we  observe  a  slight 
difference  in  the  first  few  meters  between  IRTOOL  and  IRBLEM.  The  first  derivative  of 
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refractivity  in  IRTOOL  is  bigger  than  that  in  IRBLEM.  That  means  that  the  path 
curvature  of  the  rays  in  IRTOOL  for  the  first  few  meters  is  stronger  or  in  another  words 
the  radius  of  curvature  of  the  ray  is  smaller.  So,  we  expect  that  the  MMR  (Minimum 
Mirage  Range)  of  IRTOOL  moves  toward  shorter  ranges.  This  strength  of  refraction  is 
the  reason  (for  IRTOOL)  that  the  target  at  4.32  m  is  at  the  limit  of  detection.  But  if  the 
MMR  in  IRTOOL  moves  to  shorter  ranges  than  IRBLEM,  then  we  should  expect  a 
mirage  at  least  for  the  target  at  5.77  m.  So,  we  edit  the  cut-off  height  in  IRBLEM 
(FWH=0.95)  for  the  same  reasons  that  we  mentioned  in  event  P2010F2.  The  significant 
wave  height  that  was  calculated  and  used  from  IRTOOL  was  H^O.66  m  (bigger  than  the 

measured  wave  height)  and  for  IRBLEM  was  H,=0.43  m  (smaller  than  the  measured 

wave  height).  The  results  with  the  edit  cut-off  height  of  IRBLEM  are  shown  in  Figure  34. 
Here  IRBLEM  does  not  predict  the  mirages  for  the  targets  at  5.12  m  and  5.77  m. 

Figure  35  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  08:15  to  09:40  GMT.  The  event  that  we  analyzed  took  place  from  08:08  to  08:10 
GMT  which  is  very  close  to  the  period  in  which  the  measurement  of  ASTD  were  made  by 
the  ship.  As  we  can  see  from  Figure  35  both  the  air  temperature  and  water  temperature 
increase  linearly  with  the  range.  The  maximum  variation  that  was  observed  for  ASTD 
was  about  0.3 0  C  and  so,  uncertainties  due  to  any  horizontal  inhomogeneities  can  largely 
be  neglected.  Using  also  the  fact  that  the  calculated  cut-off  height  in  IRTOOL  (0.66  m)  is 
very  close  to  the  measured  significant  wave  height  (0.5 1  m),  this  is  a  representative  event 
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in  which  it  is  shown  that  IRTOOL  calculates  a  stronger  refractive  atmosphere  than 
IRBLEM. 

4.  T2010B 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (V2)  camera  positioned  at  height  h  f  =  6. 1 1  m  and  the  infrared 

(IR1)  camera  positioned  at  height  h  t  =  6.56  m  above  the  water  level. 

In  Figure  36  the  solid  line  with  the  black  circles  represents  the  MIVR  for  a  non- 
reffactive  atmosphere  as  calculated  from  IRBLEM,  the  solid  line  with  the  white  circles 
represents  the  MIVR  for  a  non-refractive  atmosphere  as  calculated  from  IRTOOL,  the 
solid  line  with  the  squares  represents  the  observed  MIVR,  the  solid  line  with  the  rhombus 
represents  the  predicted  MIVR  from  IRBLEM  and  the  solid  line  with  the  triangles  the 
MIVR  predicted  from  the  IRTOOL  code.  The  Maximum  Intervision  Range  (MIVR) 
constitutes  for  a  given  target  height  an  absolute  limit  of  detection  range,  which  is  less 
than  the  geometrical  optical  sight  range  [Ref  14].  In  Figure  38,  the  solid  line  with  the 
squares  represents  the  observed  MMR,  the  solid  line  with  the  rhombus  represents  the 
predicted  MMR  from  IRBLEM  and  the  solid  line  with  the  triangles  the  MMR  predicted 
from  the  IRTOOL  code.  All  the  curves  that  were  used  in  Figures  36,  38  to  connect  the 
various  points  are  second  order  polynomial  fits.  The  Minimum  Mirage  Range  (MMR)  is 
the  range  that  represents  the  beginning  of  the  mirage  zone  and  is  determined  by  the  ray 
with  the  lowest  initial  elevation  which  does  not  hit  the  water  surface.  As  we  can  see  from 
Figure  36  (MIVR)  for  the  visible  V2  camera  there  is  a  difference  (about  850  m)  in  the 
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prediction  of  the  MIVR  of  a  non-refractive  atmosphere  between  the  IRTOOL  and 
IRBLEM.  This  range  corresponds  to  the  maximum  detection  range  for  a  given  target 
height  and  it  is  affected  only  by  the  earth’s  curvature  and  the  wave  height.  Equation  17 
that  gives  the  HLR  does  not  take  into  account  the  wave  height.  In  order  for  the  wave 
height  to  be  a  factor  in  the  determination  of  the  maximum  detection  range  Eq.  17  must  be 
modified  as  follows: 

HLR  =  ,/2(a  +  H.)*  (Vft  -#..)  +  ))  »  10"!  s 

V2 10->  = 

3.57  x  +  Vft-ff.  ))  (50) 

where  HLR  is  expressed  in  kilometers 

h  s  is  the  sensor  height  in  meters 

h,  is  the  target  height  in  meters 

a  is  the  earth  radius  expressed  in  meters  ( a  «  6370000  m)  and 
H  „  is  the  wave  height  in  meters. 

Using  as  Hw  in  Eq.  50  the  significant  wave  height  that  was  produced  by  IRTOOL  we 
obtain  the  same  results  as  those  calculated  from  IRTOOL  for  the  MIVR  of  a  non- 
reffactive  atmosphere  (for  the  given  target  heights).  On  the  other  hand,  the  results 
calculated  from  IRBLEM  for  the  MIVR  in  a  non-refractive  atmosphere  are  not  in 
agreement  with  Eq.  50  but  with  Eq.  17.  That  means  IRBLEM  does  not  take  into  account 
the  wave  height  (significant  wave  height)  in  the  calculation  for  the  prediction  of  the 
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geometrical  optical  sight  range.  We  also  verified  this  by  the  fact  that  when  we  edited  the 
factor  wave  height  (FWH)  using  a  much  bigger  value  than  the  nominal  the  results  for  the 
geometrical  optical  sight  range  did  not  change  at  all. 

From  Figure  36  we  observe  that  the  MIVRs  predicted  from  IRTOOL  under 
refractive  conditions  are  about  2.7  to  2.76  km  shorter  than  the  measured  values,  while  for 
IRBLEM  they  are  2.16  to  2.38  km  shorter. 

In  Figures  37A,  37B,  37C,  37D,  37E,  37F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  as  they  were 
generated  from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (symbol  “x”)  and  the  custom  atmosphere  (symbol  “o”)  and  also  those  that 
were  generated  from  the  refraction  module  of  IRBLEM  (solid  line).  From  Table  17  we 
see  that  the  air  temperature  and  humidity  were  measured  at  12  m  while  the  anemometer 
was  placed  at  a  height  of  22.6  m.  For  the  reference  height  in  IRTOOL  we  used  12  m.  This 
will  produce  an  error  less  than  2%  for  the  correct  value  of  the  wind  speed  at  the  height  of 
22.6  m,  and  so  we  do  not  expect  the  results  to  change  significantly.  But  in  order  to 
minimize  any  error  we  generated  the  wind  speed  profile  from  the  NPS  profiler.  The  value 
of  the  wind  speed  as  it  was  extracted  at  12  m  was  used  in  the  input  parameters  of 
IRTOOL.  From  the  Figures  37A  and  37B  we  observe  that  at  0  m  the  value  of  the 
humidity  is  again  greater  than  1 00%  and  the  temperature  1 0  C  less  than  the  measured 
value.  In  IRTOOL  for  the  same  height  the  relative  humidity  is  100%  and  the  water 
temperature  12.5°C  which  is  in  agreement  with  the  value  given  in  Table  17.  From  Figure 


37D  we  see  a  slight  difference  in  the  wind  speed  profiles  and  also  we  observe  that  the 
wind  speed  at  the  sea  surface  calculated  from  IRBLEM  is  0.69  m/sec.  From  Figure  37C 
we  observe  an  error  in  the  values  of  pressure  using  the  midlatitude  winter  atmosphere. 
This  error  is  not  bigger  than  0.9%  and  so  it  will  not  change  the  first  derivative  of 
refractivity.  This  is  shown  in  Figure  37F  in  which  the  two  atmospheres  that  were  used  in 
IRTOOL  give  the  same  results  for  the  gradient  of  refractivity  (ray  bending).  Also,  we 
observe  a  difference  in  the  first  few  meters  between  IRTOOL  and  IRBLEM.  The  first 
derivative  of  refractivity  in  IRTOOL  is  bigger  than  in  IRBLEM.  That  means,  the  path  of 
the  rays  in  IRTOOL  for  the  first  few  meters  is  more  curved,  or  in  other  words  the  radius 
of  curvature  of  the  ray  is  smaller.  So,  we  expect  that  the  MIVR  calculated  by  IRTOOL 
will  move  toward  shorter  ranges.  This  is  the  reason  for  the  differences  between  the 
results. 

From  Figure  38  we  observe  that  IRTOOL  calculates  a  MMR  for  the  given  target 
heights  which  is  about  1.44  to  2.5  km  longer  than  the  observed  one,  while  IRTOOL 
calculates  a  MMR  about  0.4  to  0.6  km  shorter.  The  significant  wave  height  that  was  used 
is  a  very  important  (crucial)  factor  for  the  determination  of  the  MMR.  We  expect  that  an 
increased  significant  wave  height  would  move  the  predicted  MMRs  to  longer  ranges  for 
both  models  and  also  would  lower  the  maximum  predicted  height  for  a  mirage.  The 
significant  wave  height  calculated  from  IRTOOL  was  0.61  m  while  for  IRBLEM  it  was 
0.1626  m.  So,  in  order  to  compare  both  models  for  the  same  cut-off  wave  height  we 
edited  the  Factor  Wave  Height  (FWH)  in  IRBLEM.  The  reason  that  we  edited  the  height 
in  IRBLEM  instead  in  IRTOOL  has  been  explained  in  the  analysis  of  the  event  P2010F2. 
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The  MMR  calculated  from  IRBLEM  which  is  now  about  1. 8-3.2  km  longer  than  the 
recorded  one  is  shown  in  Figure  38  (solid  line  with  the  symbol  The  change  in  the 
wave  height  did  not  change  the  results  for  the  MIVR  in  IRBLEM.  This  is  because  the 
MIVR  is  not  affected  as  much  as  the  MMR  from  the  change  of  the  significant  wave 
height  unless  we  use  a  very  big  wave  height. 

Figure  25  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  during  that  event. 
The  ASTD  that  we  used  was  -4.1  °C  but  according  to  Figure  25  there  is  a  variation  in 
ASTD  of  approximately  2°C  from  10-20  km.  That  means  the  results  can  be  further 
improved  if  we  take  into  consideration  this  figure.  Using  a  smaller  magnitude  of  ASTD 
the  strength  of  the  refraction  for  both  models  will  be  reduced,  resulting  in  curves  that  will 
move  to  greater  ranges  i.e.,  toward  the  measured  data  (under  the  assumption  that  both 
models  use  the  same  significant  wave  height  of  0.1626  m). 

Figures  39,  41  show  the  results  of  the  model  comparison  for  the  data  observed  by 
the  IR1  (3-5  //m)  infrared  camera  positioned  at  height  ht  =  6.56  m  above  the  water 

level.  As  we  can  see  from  Figure  39  (MIVR)  the  IRTOOL  model  results  for  the  MIVR 
are  2.3  to  2.5  km  shorter  and  the  IRBLEM  results  are  1.6  to  2.1  km  shorter  than  the 
observed  ones.  Also  the  non-refraction  range  that  is  predicted  from  IRBLEM  is  0.7- 1.0 
km  longer  than  the  one  predicted  from  IRTOOL.  That  happens  because  IRBLEM  did  not 
take  into  account  the  wave  height  for  the  prediction  of  the  geometrical  optical  sight  range 
as  we  have  already  mentioned  in  the  analysis  of  this  event  with  the  visible  (V2)  camera. 
The  non-refraction  range  predicted  from  IRTOOL  is  shorter  than  the  real  one  (if  we  could 
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have  measured  it)  because  the  significant  wave  height  that  was  used  in  Eq.  50  (0.61  m)  is 
bigger  than  the  real  measured  one  (0.17  m).  All  the  profiles  calculated  from  IRTOOL  are 
in  agreement  with  the  profiles  that  were  calculated  for  the  same  event  with  the  V2  camera 
because  IRTOOL  does  not  model  dispersion.  For  IRBLEM  the  profiles  of  temperature, 
humidity,  pressure  and  wind  speed  remain  the  same  while  the  profiles  of  refractivity 
(Figure  40A)  and  its  gradient  (Figure  40B)  differ  slightly  from  before,  due  to  the  different 
coefficients  A,  B  that  are  used  for  the  mid-infrared  in  the  equation  of  refractivity.  So,  we 
expect  about  the  same  results  as  before.  From  Figure  41  we  observe  that  the  IRTOOL 
model  results  for  the  MMR  are  1 .65  to  2. 1 7  km  longer  and  the  IRBLEM  results  are  0.4  to 
0.8  km  shorter  than  the  observed  ones.  The  significant  wave  height  that  was  used  for  the 
prediction  of  those  results  is  0.61  m  for  IRTOOL  and  0.1626  m  for  IRBLEM.  When  we 
edited  the  FWH  of  IRBLEM  to  produce  the  same  wave  height  as  IRTOOL  the  MMR 
calculated  from  IRBLEM  was  1 .94  to  2.7  km  toward  greater  ranges  than  the  observed 
values.  The  results  are  shown  in  Figure  42  (solid  line  with  the  symbol  Also,  this 
change  in  the  wave  height  did  not  alter  the  results  for  the  MIVR  produced  from  IRBLEM 
at  all  for  the  reason  that  we  explained  earlier. 

Figure  25  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  during  that  event. 
The  ASTD  that  we  used  was  -4.1  °C  but  according  to  Figure  25  there  is  a  variation  in 
ASTD  of  approximately  2°C  from  10-20  km.  That  means  the  results  can  be  further 
improved  if  we  take  into  consideration  this  figure.  Using  a  smaller  magnitude  of  ASTD 
the  strength  of  the  refraction  for  both  models  will  be  reduced,  resulting  in  curves  that  will 
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move  to  greater  ranges  i.e.,  toward  the  measured  data  (under  the  assumption  that  both 
models  use  the  same  significant  wave  height  of  0.16  m).  Taking  into  account  the  fact  that 
the  equation  used  for  the  refractivity  in  IRTOOL  and  the  equation  used  in  IRBLEM  for 
the  mid-infrared  are  almost  identical  and  also  the  fact  that  the  pressure  and  humidity 
profiles  are  in  good  agreement  we  can  say  that  the  difference  in  the  results  occurs  due  to 
the  difference  in  the  temperature  profiles  between  the  two  models  for  the  first  few  meters. 

The  use  of  the  custom  atmosphere  profile  in  IRTOOL  instead  of  the  mid-latitude 
winter  atmosphere  does  not  change  the  results  at  all  as  we  saw  for  all  the  events  that  we 
examined  until  now.  For  that  reason,  for  the  analysis  of  the  following  events  the  mid¬ 
latitude  winter  atmosphere  will  be  used. 

5.  P2210J2 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (V2)  positioned  at  height  hs  =  8.421  m  above  the  water 
level. 

From  Table  17  we  see  that  all  the  measurements  were  made  at  the  same  height  of 
3.4  m.  As  we  can  see  from  the  observed  data  in  Figure  42,  a  mirage  was  observed  for  the 
targets  at  3.82  and  4.52  m  but  not  for  the  target  at  3.07  m.  Both  models  IRBLEM  and 
IRTOOL  do  not  detect  the  target  at  3.07  m,  and  also  they  do  not  predict  mirage  for  the 
targets  at  3.82  m  and  4.52  m.  For  a  non-reffactive  atmosphere  any  target  below  «2  m 
should  be  below  the  horizon.  The  results  that  both  models  give  are  in  good  agreement 
with  each  other  but  not  with  the  observed  data. 
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In  Figures  43A,  43B,  43C,  43D,  43E  and  43F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  reffactivity  and  gradient  of  refractivity  as  they  were 
generated  from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (dash  line)  and  those  that  were  generated  from  the  refraction  module 
(LWKD)  of  IRBLEM  (solid  line).  From  the  Figures  43A  and  43B  we  observe  that  at  0  m 
the  value  of  the  humidity  is  again  greater  than  100%  and  the  temperature  1.4  °C  less  than 
the  measured  one.  In  IRTOOL  for  the  same  height  the  relative  humidity  is  100%  and  the 
water  temperature  12  °C  which  is  in  agreement  with  the  value  given  in  Table  17.  From 
Figure  43D  we  observe  very  good  agreement  in  the  wind  speed  profiles  and  also  that  the 
wind  speed  at  the  water  surface  calculated  from  IRBLEM  is  now  0.01  m/sec.  From  Figure 
43  C  the  error  in  the  pressure  profile  using  the  midlatitude  winter  atmosphere  is  not 
greater  than  0.97%,  and  so  we  do  not  expect  any  change  in  the  results.  Because  neither 
model  predicts  a  mirage  for  any  target  we  expect  that  the  results  will  not  change 
significantly  when  we  edit  the  cut-off  height  in  IRBLEM.  The  significant  wave  height 
that  was  calculated  and  used  from  IRTOOL  was  Hs=2.29  m  (much  bigger  than  the  real 

measured  wave  height)  and  for  IRBLEM  was  H^O.76  m  (smaller  than  the  real  measured 

wave  height).  The  results  with  the  edit  cut-off  height  of  IRBLEM  are  shown  in  Figure  44 
and  as  we  see  they  do  not  change  significantly  (the  changes  can  not  be  observed  on  the 
scale  we  use). 

Figure  45  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  20:00  to  22:10  GMT.  The  event  that  we  analyzed  took  place  from  20:22  to  20:24 
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GMT  which  is  close  to  the  period  in  which  the  measurements  of  ASTD  were  made  by  the 
ship.  As  we  can  see  from  Figure  45  there  is  a  variation  in  ASTD  of  about  0.8 0  C  and  so 
this  horizontal  inhomogeneity  could  have  an  effect  on  the  agreement  between  the  models 
and  the  data. 

6.  P2510M1 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (VI)  positioned  at  height  h  5  =  7.51  m  above  the  water  level 

and  also  the  visible  (V2)  camera  positioned  at  height  h  s  =  4.32  m  above  the  water  level. 

From  Table  17  we  see  that  the  all  the  meteorological  measurements  were  made  at 
the  same  height  of  3.4  m.  As  we  can  see  from  the  recorded  data  in  Figure  46,  no  mirage 
was  observed  for  the  given  target  heights  and  both  models  IRBLEM  and  IRTOOL  do  not 
predict  any  mirage  for  any  target  height  for  this  set  of  conditions.  For  a  non-reffactive 
atmosphere  any  target  below  «2  m  should  be  below  the  horizon.  The  results  that  both 
models  give  are  in  good  agreement  among  themselves  and  also  with  the  experimental 
data.  In  Figures  47A,  47B,  47C,  47D,  47E  and  47F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  as  generated 
from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (dash  line)  and  those  that  were  generated  from  the  refraction  module 
(LWKD)  of  IRBLEM  (solid  line).  From  the  Figures  47A  and  47B  we  observe  that  at  0  m 
the  value  of  the  humidity  is  now  about  100%  (100.92%)  and  the  temperature  11.3°C 
which  is  almost  in  agreement  with  the  measured  value.  In  IRTOOL  for  the  same  height 
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the  relative  humidity  is  100%  and  the  water  temperature  11.6°C  in  agreement  with  the 
value  given  in  Table  17.  From  Figure  47D  we  observe  very  good  agreement  in  the  wind 
speed  profiles  and  also  that  the  wind  speed  at  the  water  surface  calculated  from  IRBLEM 
is  now  0.01  m/sec.  Looking  in  the  Figure  47C  the  error  in  the  pressure  profile  using  the 
midlatitude  winter  atmosphere  is  not  greater  than  0.97%,  and  so  we  do  not  expect  any 
change  in  the  results.  In  Figure  47E  we  observe  that  the  gradient  of  refractivity  is  almost 
identical  between  the  two  models  and  also  the  temperature  profiles  (Figure  47B),  with  a 
slight  difference  at  heights  above  5  m.  Because  neither  model  predicts  a  mirage  for  any 
target  we  expect  that  the  results  are  not  going  to  change  significantly  when  we  edit  the 
cut-off  height  in  IRBLEM.  The  significant  wave  height  that  was  calculated  and  used  from 
IRTOOL  was  =1.75  m  (much  bigger  than  the  real  measured  wave  height)  and  for 

IRBLEM  was  H  s  =0.73  m  (smaller  than  the  real  measured  wave  height).  The  results  with 
the  edit  cut-off  height  of  IRBLEM  (not  shown)  are  the  same  as  in  Figure  46. 

Figure  48  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  03:10  to  05:25  GMT.  The  event  that  we  analyzed  took  place  from  03:17  to  03:18 
GMT.  As  we  can  see  from  Figure  48  there  is  a  variation  in  ASTD  of  about  0.8 0  C,  and  so 
this  horizontal  inhomogeneity  could  have  an  effect  on  the  agreement  between  the  models 
and  the  data. 

With  the  visible  (V2)  camera  the  results  are  the  same  as  with  the  visible  (VI) 
camera.  From  Table  17  we  see  that  all  the  measurements  were  made  at  the  same  height  of 
3.4  m.  As  we  can  see  in  Figure  49  no  mirage  was  observed  for  the  given  target  heights 
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and  neither  model  IRBLEM  or  IRTOOL  predicts  any  mirage  for  any  target  height  for  this 
set  of  conditions.  For  a  non-refractive  atmosphere  any  target  below  «3.5  m  should  be 
below  the  horizon.  The  results  that  the  two  models  give  are  in  good  agreement  between 
themselves  and  also  with  the  experimental  data.  The  profiles  of  humidity,  temperature, 
pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  for  IRTOOL  and  IRBLEM 
are  the  same  as  the  previous  ones  (47A-47F),  because  IRTOOL  does  not  model 
dispersion  and  because  we  are  still  in  the  visible,  so  the  coefficients  A,  B  of  IRBLEM  do 
not  change.  The  significant  wave  heights  that  were  used  from  both  models  are  the  same 
as  previously  and  the  results  when  we  edit  the  cut-off  height  in  IRBLEM  will  not  change. 

7.  P2810P3 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
camera  used  was  the  visible  (V2)  positioned  at  height  hv  =  4.64  m  above  the  water  level. 

From  Table  17  we  see  that  the  all  the  measurements  were  made  at  the  same  height 
of  3.4  m.  As  we  can  see  from  Figure  50,  neither  model  IRBLEM  or  IRTOOL  predicts  the 
mirage  for  the  target  at  5.81  m  but  the  rest  of  the  results  are  in  very  good  agreement 
between  themselves  and  with  the  recorded  data.  For  a  non-refractive  atmosphere  any 
target  below  « 1  m  should  be  below  the  horizon.  This  height  increases  to  about  ~4  m  for 
IRBLEM  and  a  little  more  for  IRTOOL  (4.2  m).  This  is  a  very  representative  event 
because  the  cut-off  height  calculated  from  IRTOOL  was  0.49  m,  almost  the  same  as  the 
real  measured  significant  wave  height  (0.51  m).  Also,  the  cut-off  height  that  was  used 
from  IRBLEM  was  0.43  m.  That  means,  for  this  event  both  models  were  compared  under 
the  same  (almost)  cut-off  height,  very  close  to  the  real  one,  and  so  the  results  that  we  get 
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are  representative  of  their  ability  to  predict  the  values  of  the  real  data.  Even  if  we  edit  the 
cut-off  height  in  IRBLEM,  so  as  to  agree  with  the  cut-off  height  of  IRTOOL,  the  results 
in  IRBLEM  will  not  change  (any  change  will  not  be  observable  on  the  range  scale  that  we 
use)  because  of  the  small  difference  in  the  cut-off  heights  between  the  models. 

In  Figures  51A,  51B,  51C,  51D,  51E  and  51F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  as  generated 
from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (dash  line)  and  those  that  were  generated  from  the  refraction  module 
(LWKD)  of  IRBLEM  (solid  line).  From  Figures  51A  and  51B  we  observe  that  at  0  m  the 
value  of  the  humidity  is  104.8%  which  is  not  a  physical  result,  and  the  temperature  0.8 0  C 
less  than  the  measured  value,  in  spite  of  the  fact  that  we  used  the  measured  value  in  the 
inputs  to  IRBLEM.  In  IRTOOL  for  the  same  height  the  relative  humidity  is  100%  and  the 
water  temperature  11.7°C  in  agreement  with  the  value  given  in  Table  17  and  the  value 
that  we  used  in  the  IRTOOL  inputs.  From  Figure  5 ID  we  observe  very  good  agreement  in 
the  wind  speed  profiles  and  also  the  wind  speed  at  the  water  surface  calculated  from 
IRBLEM  is  0.72  m/sec.  Looking  in  Figure  51C  the  error  in  the  pressure  profile  using  the 
midlatitude  winter  atmosphere  is  not  greater  than  0.91%,  and  so  we  do  not  expect  any 
change  in  the  results.  In  Figure  51 F  we  observe  that  the  gradients  of  refractivity  for  the 
two  models  are  in  very  good  agreement.  For  that  reason  both  models  predict  almost  the 
same  results. 
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Figure  35  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  08:15  to  09:40  GMT.  The  event  that  we  analyzed  took  place  from  09:18  to  09:24 
GMT.  As  we  can  see  from  Figure  35  both  the  air  temperature  and  water  temperature 
increase  linearly  with  the  range.  The  maximum  variation  that  we  measured  for  ASTD  was 
about  0.3  0  C,  and  so  uncertainties  due  to  any  horizontal  inhomogeneities  can  largely  be 
neglected.  As  a  result  this  is  also  a  very  good  event  for  the  comparison  between  the  two 
models. 

8.  T2810K 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
cameras  used  were  the  visible  (VI)  camera  and  the  infrared  (IR1)  camera  both  positioned 
at  height  h,  =  7.83  m  above  the  water  level. 

In  Figure  52  the  black  solid  line  with  the  black  circles  represents  the  MIVR  for  a 
non-reif active  atmosphere  as  calculated  from  IRJBLEM,  the  black  solid  line  with  the 
white  circles  represents  the  MIVR  for  a  non-refractive  atmosphere  as  calculated  from 
IRTOOL,  the  black  solid  line  with  the  squares  represents  the  observed  MIVR,  the  black 
solid  line  with  the  rhombus  represents  the  predicted  MIVR  from  IRBLEM  and  the  black 
solid  line  with  the  triangles  the  MIVR  predicted  from  the  IRTOOL  code.  In  Figure  54,  the 
squares  represent  the  observed  MMR,  the  rhombus  represents  the  MMR  predicted  from 
IRBLEM  and  the  triangles  the  MMR  predicted  from  the  IRTOOL  code.  All  the  curves 
that  were  used  in  Figure  52  to  connect  the  various  points  are  second  order  polynomial 
fits.  As  we  can  see  from  Figure  52  (MIVR)  for  the  visible  VI  camera  there  is  a  difference 
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(about  700  m)  in  the  prediction  of  the  MIVR  in  a  non-refractive  atmosphere  between  the 
IRTOOL  and  IRBLEM.  Using  as  Hw  in  Eq.  50  the  significant  wave  height  that  was 
produced  by  IRTOOL  we  obtain  the  same  results  as  those  calculated  from  IRTOOL  for 
the  MIVR  of  a  non-refractive  atmosphere  (for  the  given  target  heights).  On  the  other 
hand,  the  results  calculated  from  IRBLEM  for  the  MIVR  in  a  non-refractive  atmosphere 
are  not  in  agreement  with  Eq.  50  but  with  Eq.  17.  That  means  IRBLEM  does  not  take  into 
account  the  wave  height  (significant  wave  height)  in  the  calculation  for  the  prediction  of 
the  geometrical  optical  sight  range. 

From  Figure  52  we  observe  that  the  MIVR  predicted  from  IRTOOL  under 
refractive  conditions  is  about  1 .26  to  1 .69  km  shorter  than  the  measured  one,  while  for 
IRBLEM  it  is  0.91  to  1 .53  km  shorter  than  measured. 

In  Figures  53A,  53B,  53C,  53D,  53E,  53F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  as  they  were 
generated  from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (dash  line)  and  also  those  that  were  generated  from  the  refraction  module  of 
IRBLEM  (solid  line).  From  the  Figures  53A  and  53B  we  observe  that  at  0  m  the  value  of 
the  humidity  is  again  greater  than  100%  and  the  temperature  0.7  °C  less  than  the 
measured  one.  In  IRTOOL  for  the  same  height  the  relative  humidity  is  100%  and  the 
water  temperature  1 1.7  °C  which  is  in  agreement  with  the  value  given  in  Table  17.  From 
Figure  53D  we  see  a  slight  difference  in  the  wind  speed  profiles  and  also  we  observe  that 
the  wind  speed  at  the  sea  surface  calculated  from  IRBLEB  is  0.72  m/sec.  From  Figure 
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53C  we  observe  an  error  in  the  values  of  pressure  using  the  midlatitude  winter 
atmosphere.  This  error  is  not  bigger  than  0.9%  and  so  it  will  not  change  the  first 
derivative  of  refractivity.  In  Figure  53F  we  observe  that  the  gradients  of  refractivity  are  in 
very  good  agreement  between  the  two  models,  although  there  is  a  slight  difference  for  the 
first  few  meters  between  the  two  models.  The  first  derivative  of  refractivity  in  IRTOOL  is 
bigger  than  in  IRBLEM.  That  means,  the  path  curvature  of  the  rays  in  IRTOOL  for  the 
first  few  meters  is  greater  or  in  other  words  the  radius  of  curvature  of  the  ray  is  smaller. 
So,  we  expect  that  the  MIVR  calculated  from  IRTOOL  will  move  toward  shorter  ranges. 
This  is  the  reason  for  the  difference  between  the  results. 

From  Figure  54  we  observe  that  IRTOOL  calculates  a  MMR  for  the  given  target 
heights  which  is  about  38  to  204  m  longer  than  the  observed  one,  while  IRBLEM 
calculates  a  MMR  about  19  to  161m  longer  than  observed.  The  significant  wave  height 
that  was  used  is  a  very  important  factor  for  the  determination  of  the  MMR.  We  expect 
that  an  increased  significant  wave  height  would  move  the  predicted  MMRs  to  longer 
ranges  for  both  models  and  also  would  lower  the  maximum  predicted  height  for  a  mirage. 
The  significant  wave  height  calculated  from  IRTOOL  was  0.49  m  while  for  IRBLEM  it 
was  0.42  m. 

This  is  an  extremely  good  case  because  both  models  use  about  the  same  cut-off 
height  which  is  also  very  close  to  the  measured  one  (0.51  m).  And  in  this  case  we  got 
very  good  results  for  both  models.  In  order  to  see  how  the  results  in  IRBLEM  will 
change  when  we  change  its  cut-off  height  we  edited  the  Factor  Wave  Height  (FWH)  so  as 
to  produce  the  same  cut-off  height  as  IRTOOL.  The  MMR  calculated  now  from 
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IRBLEM,  which  is  about  0.23-0.5  km  longer  than  .the  recorded  one,  is  shown  in  Figure  54 
(symbol  “4c”).  The  change  in  the  wave  height  did  not  change  the  results  for  the  MIVR  in 
IRBLEM.  This  is  because  the  MIVR  is  not  affected  as  much  as  the  MMR  from  the 
change  of  the  significant  wave  height  unless  a  very  big  wave  height  is  used. 

Figure  35  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  during  that 
event.  The  ASTD  that  we  used  was  -2.5  °C.  As  we  can  see  from  Figure  36  both  the  air 
temperature  and  water  temperature  increase  linearly  with  the  range.  The  maximum 
variation  that  we  measured  for  ASTD  was  about  0.3  0  C,  and  so  uncertainties  due  to  any 
horizontal  inhomogeneities  can  largely  be  neglected.  As  a  result  this  is  a  very  good  event 
for  the  comparison  between  the  two  models. 

For  the  IR1  camera  IRTOOL  predicts  exactly  the  same  results  as  with  the  visible 
(VI)  camera  because  both  cameras  were  placed  at  the  same  height  and  because  IRTOOL 
does  not  model  dispersion.  The  results  for  IRBLEM  we  expect  to  be  a  little  different 
because  the  coefficients  A,  B  that  are  used  for  the  calculation  of  refractivity  are  band 
dependent.  All  the  profiles  for  IRTOOL  remain  the  same  as  before  while  for  IRBLEM 
only  the  profiles  of  the  refractivity  and  its  first  gradient  change  (Figures  56A  and  56B). 
Figures  55,  57  show  the  results  of  the  model  comparison  for  the  MMR  and  MIVR.  As  we 
can  see  in  Figure  55  (MIVR)  the  IRBLEM  results  are  0.63  to  1.44  km  shorter  than  the 
observed  values  while  for  IRTOOL  are  0.98-1.58  km  shorter.  From  Figure  57  we  observe 
that  the  IRTOOL  and  IRBLEM  model  results  for  the  MMR  are  in  very  good  agreement 
with  the  observed  ones.  In  order  to  see  how  the  results  in  IRBLEM  will  change  when  we 
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change  its  cut-off  height  we  edited  the  Factor  Wave  Height  (FWH)  so  as  to  produce  the 
same  cut-off  height  as  IRTOOL.  The  MMR  calculated  now  from  IRBLEM,  which  is 
about  0.21-0.4  km  longer  than  the  recorded  one,  is  shown  in  Figure  57  (symbol  The 
change  in  the  wave  height  did  not  change  the  results  for  the  MIYR  in  IRBLEM. 

Figure  35  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  during  that 
event.  The  ASTD  that  we  used  was  -2.5  °C.  As  we  can  see  from  Figure  35  both  the  air 
temperature  and  water  temperature  increase  linearly  with  the  range.  The  maximum 

variation  that  we  measured  for  ASTD  was  about  0.3 0  C,  and  so  uncertainties  due  to  any 
horizontal  inhomogeneity  can  largely  be  neglected.  As  a  result  this  is  a  very  good  event 
for  the  comparison  between  the  two  models. 

9.  T2510F 

The  conditions  under  which  this  event  took  place  are  shown  in  Table  17.  The 
cameras  used  were  the  visible  (VI),  the  visible  (V2),  the  visible  (V3),  the  visible  (V4), 
the  infrared  (IR1)  and  the  infrared  (IR2)  camera.  For  all  the  plots,  the  black  solid  line  with 
the  black  circles  represents  the  MIVR  for  a  non-refractive  atmosphere  as  calculated  from 
IRBLEM,  the  black  solid  line  with  the  white  circles  represents  the  MIVR  for  a  non- 
refractive  atmosphere  as  calculated  from  IRTOOL,  the  squares  represent  the  observed 
MIVR  or  MMR,  the  rhombus  represents  the  MIVR  or  MMR  predicted  from  IRBLEM 
and  the  triangles  the  MIVR  or  MMR  predicted  from  the  IRTOOL  code. 

The  visible  (VI)  camera  was  positioned  at  h  v  =  7.55  m  above  the  water  level.  As 
we  can  see  from  Figure  58  there  is  a  difference  (about  2  km)  in  the  prediction  of  the 
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MIVR  in  a  non-refractive  atmosphere  between  the  IRTOOL  and  IRBLEM.  Using  as  HM, 

in  Eq.  50  the  significant  wave  height  that  was  produced  by  IRTOOL  (1.46  m)  we  obtain 
the  same  results  as  those  calculated  from  IRTOOL  for  the  MIVR  of  a  non-refractive 
atmosphere  (for  the  given  target  heights).  On  the  other  hand,  the  results  calculated  from 
IRBLEM  for  the  MIVR  in  a  non-refractive  atmosphere  are  not  in  agreement  with  Eq.  50 
but  with  Eq.  17.  That  means  IRBLEM  does  not  take  into  account  the  wave  height 
(significant  wave  height)  in  the  calculation  for  the  prediction  of  the  geometrical  optical 
sight  range.  From  the  same  figure  we  observe  that  the  MIVR  predicted  from  IRTOOL 
under  refractive  conditions  is  about  0.85  to  1.65  km  shorter  than  the  measured  value, 
while  for  IRBLEM  it  is  700  to  950  m  shorter.  It  is  interesting  to  note  that  even  though  we 
have  subrefractive  conditions  (ASTD  =  -0.6  °C)  the  MIVR  is  at  greater  ranges  than  the 
geometrical  horizon.  The  ray  curvature  depends  upon  the  index  of  refraction  gradient, 
which  depends  upon  the  gradients  in  density.  For  neutral  stability  conditions  (ASTD 
about  0),  the  air  density  is  still  decreasing  with  height  and  the  ray  is  curved  around  the 
earth  slightly  compared  to  a  straight  ray.  It  apparently  takes  a  more  negative  ASTD  than  - 
0.6  °C  (depending  also  on  wind  speed)  to  undo  this  effect. 

From  the  profile  of  the  air  density  (not  shown)  we  observe  that  the  gradient  of  the 
air  density  became  negative  at  the  height  of  about  2  m.  That  means  the  air  density  starts  to 
decrease  at  that  height.  Taking  into  account  that  the  significant  wave  height  produced  by 
IRTOOL  was  1.46  m  and  also  the  fact  that  the  sensor  was  placed  at  7.55  m  above  the 
water  level  we  have  a  height  region  in  which  the  air  density  decreases  with  height.  So,  the 
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rays  bend  downward  and  we  detect  at  ranges  beyond  the  horizon.  This  can  also  be  seen 
from  the  gradient  of  refractivity  (Figure  59F)  in  which  the  gradient  starts  to  become 
negative  at  the  height  of  2  m. 

In  Figures  59A,  59B,  59C,  59D,  59E,  59F  we  see  the  profiles  of  humidity, 
temperature,  pressure,  wind  speed,  refractivity  and  gradient  of  refractivity  as  generated 
from  the  atmosphere  profiler  module  of  IRTOOL  using  the  midlatitude  winter 
atmosphere  (dash  line)  and  also  those  that  were  generated  from  the  refraction  module  of 
IRBLEM  (solid  line).  From  Figures  59A  and  59B  we  observe  that  at  0  m  the  value  of  the 
humidity  is  102%  which  is  a  non-physical  result,  and  the  temperature  0.24  °C  less  than 
the  measured  value.  In  IRTOOL  for  the  same  height  the  relative  humidity  is  100%  (sea 
surface  100%  saturated)  and  the  water  temperature  11.6°C  which  is  in  agreement  with 
the  value  given  in  Table  17.  From  Figure  59D  we  see  a  slight  difference  in  the  wind 
speed  profiles  above  5  m.  From  Figure  59C  we  observe  an  error  in  the  values  of  pressure 
using  the  midlatitude  winter  atmosphere.  This  error  is  not  bigger  than  1 .2%  and  so  it  will 
not  alter  significantly  the  first  derivative  of  refractivity.  In  Figure  59F  we  observe  that  the 
gradients  of  refractivity  are  in  very  good  agreement  between  the  two  models,  although 
there  is  a  slight  difference  for  the  first  few  meters  between  the  two  models.  The  first 
derivative  of  refractivity  in  IRTOOL  is  bigger  than  in  IRBLEM.  That  means,  the  path 
curvature  of  the  rays  in  IRTOOL  for  the  first  few  meters  is  greater,  or  in  other  words  the 
radius  of  curvature  of  the  ray  is  smaller.  So,  we  expect  that  the  MIYR  calculated  with 
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IRTOOL  will  move  toward  shorter  ranges.  The  results  that  are  shown  in  Figure  58  are  in 
very  good  agreement  between  both  models  and  the  measured  data. 

In  Figure  60  we  observe  that  MMRs  were  observed  for  targets  from  1 1  m  to  21  m 
above  the  water  level  but  not  at  3.34  m,  6.08m  and  6.38  m  above  the  water  level. 
IRTOOL  does  not  predict  a  mirage  for  the  given  target  heights  while  IRBLEM  predicts 
the  existence  of  weak  secondary  images,  or  image  stretching,  for  all  the  targets  above  6 
m.  The  relative  strengths  for  the  formation  of  secondary  images  can  be  seen  by  comparing 
the  range  difference,  at  any  elevation,  between  the  predicted  MIVR  and  MMR. 

The  significant  wave  height  that  was  used  is  a  very  important  factor  for  the 
determination  of  the  MMR.  We  expect  that  an  increased  significant  wave  height  would 
move  the  predicted  MMRs  to  longer  ranges  for  both  models  and  also  would  lower  the 
maximum  predicted  height  for  a  mirage.  The  significant  wave  height  calculated  from 
IRTOOL  was  1.46  m  while  for  IRBLEM  it  was  0.74  m.  In  order  to  see  how  the  results  in 
IRBLEM  will  change  when  we  change  its  cut-off  height  we  edited  the  Factor  Wave 
Height  (FWH)  so  as  to  produce  the  same  cut-off  height  as  IRTOOL.  Now  IRBLEM  does 
not  predict  mirages  for  the  given  target  heights.  The  change  in  the  wave  height  changed 
slightly  the  results  for  the  MIVR  in  IRBLEM.  Now  the  MIVRs  move  towards  shorter 
ranges  (Figure  58,  symbol 

The  visible  (V2)  camera  was  positioned  at  hs  =  9.34  m  above  the  water  level.  As 

we  can  see  from  Figure  61  the  MIVR  predicted  from  IRTOOL  under  refractive  conditions 
is  about  1  to  1.4  km  shorter  than  the  measured  value,  while  for  IRBLEM  it  is  0.7  to  1  km 
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shorter.  The  profiles  of  humidity,  temperature,  pressure,  wind  speed,  refractivity  and 
gradient  of  refractivity  for  both  models  are  identical  to  those  that  were  generated  for  the 
same  event  with  the  visible  (VI)  camera  (59A-59F).  In  Figure  62  we  observe  that  MMRs 
were  measured  only  for  the  target  at  20.7  m  above  the  water  level.  IRTOOL  does  not 
predict  a  mirage  for  the  given  target  heights  while  IRBLEM  predicts  the  existence  of 
weak  secondary  images  for  all  targets  above  6  m.  The  mirage  range  for  the  target  at  20.7 
m  is  about  750  m  shorter  than  the  measured  value.  The  significant  wave  height  that  was 
used  is  a  very  important  factor  for  the  determination  of  the  MMR.  We  expect  that  an 
increased  significant  wave  height  would  move  the  predicted  MMRs  to  longer  ranges  for 
both  models  and  also  would  lower  the  maximum  predicted  height  for  a  mirage.  The 
significant  wave  height  calculated  from  IRTOOL  was  1 .46  m  while  for  IRBLEM  it  was 
0.74  m.  In  order  to  see  how  the  results  in  IRBLEM  will  change  when  we  change  its  cut¬ 
off  height  we  edited  the  Factor  Wave  Height  (FWH)  so  as  to  produce  the  same  cut-off 
height  as  IRTOOL.  Now,  neither  IRBLEM  nor  IRTOOL  predicts  a  mirage  for  the  given 
target  heights.  The  change  in  the  wave  height  changed  slightly  the  results  for  the  MIVR  in 
IRBLEM  (Figure  61,  symbol 

The  visible  (V3)  camera  was  positioned  at  h  v  =  20.51  m  above  the  water  level.  As 

we  can  see  from  Figure  63  the  MIVR  predicted  from  IRTOOL  under  refractive  conditions 
is  about  2.3  to  2.9  km  longer  than  measured,  while  for  IRBLEM  it  is  2.88  to  3  km  longer. 
The  profiles  of  humidity,  temperature,  pressure,  wind  speed,  refractivity  and  gradient  of 
refractivity  for  both  models  are  the  same  as  those  that  were  generated  for  the  same  event 
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with  the  visible  (VI)  camera  (59A-59F).  Figure  64  shows  that  IRTOOL  does  not  predict 
the  existence  of  any  mirages  and  also  that  no  MMR  was  observed  for  the  given  target 
heights,  while  IRBLEM  predicts  the  existence  of  weak  secondary  images  for  all  targets 
above  1 1  m.  The  significant  wave  height  that  was  used  is  a  very  important  factor  for  the 
determination  of  the  MMR.  The  significant  wave  height  calculated  from  IRTOOL  was 
1.46  m  while  for  IRBLEM  it  was  0.74  m.  In  order  to  see  how  the  results  in  IRBLEM  will 
change  when  we  change  its  cut-off  height  we  edited  the  Factor  Wave  Height  (FWH)  so  as 
to  produce  the  same  cut-off  height  as  IRTOOL.  Now  IRBLEM  does  not  predict  a  mirage 
for  the  given  target  heights.  The  change  in  the  wave  height  changed  slightly  the  results  for 
the  MIVR  in  IRBLEM  (Figure  63,  symbol 

The  visible  (V4)  camera  was  positioned  at  hf=  7.14  m  above  the  water  level.  As 

we  can  see  from  Figure  65  the  MIVR  predicted  from  IRTOOL  under  refractive  conditions 
is  about  0.5  to  1 .4  km  shorter  than  the  measured  one,  while  for  IRBLEM  it  is  0.4  to  0.65 
km  shorter.  The  profiles  of  humidity,  temperature,  pressure,  wind  speed,  refractivity  and 
gradient  of  refractivity  for  both  models  are  the  same  as  those  that  were  generated  for  the 
same  event  with  the  visible  (VI)  camera  (59A-59F).  Figure  66  shows  that  only  a  single 
MMR  data  point  was  observed,  due  to  the  target  at  height  20.7  m,  and  also  that  neither 
model  predicted  the  existence  of  any  secondaiy  images.  When  we  edited  the  cut-off 
height  in  IRBLEM  so  as  to  be  the  same  as  in  IRTOOL  the  results  related  to  MMR  did  not 
change  as  we  expected,  but  the  MIVRs  were  moved  a  little  (Figure  65,  symbol  “*”). 
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The  infrared  (IR1)  camera  was  positioned  at  hi.=  7.55  m  above  the  water  level.  As 


we  can  see  from  Figure  67  the  MIVR  predicted  from  IRTOOL  under  refractive  conditions 
is  about  0.7  to  1.5  km  shorter  than  the  measured  one,  while  for  IRBLEM  it  is  0.5  to  0.8 
km  shorter.  All  the  profiles  for  IRTOOL  remain  the  same  as  before  while  for  IRBLEM 
only  the  profiles  of  the  refractivity  and  its  first  gradient  changed  (Figures  68A  and  68B) 
due  to  the  band  dependent  coefficients  A,  B  that  are  used  for  the  calculation  of 
refractivity.  As  we  can  see  from  Figure  68B  the  curves  of  the  gradient  of  refractivity  for 
the  two  .models  are  in  better  agreement  compared  with  the  visible  cameras.  That 
happened  because  in  the  mid-infrared  the  equations  that  are  used  in  the  two  models  for 
the  calculation  of  refractivity  are  approximately  the  same.  Figure  69  shows  that  IRTOOL 
does  not  predict  the  existence  of  any  mirages  for  the  given  target  heights  and  also  that 
only  one  MMR  data  point  was  observed,  due  to  the  target  at  20.7  m.  IRBLEM  predicts 
the  existence  of  weak  secondary  images  for  all  targets  above  6  m.  The  significant  wave 
height  that  was  used  is  a  very  important  factor  for  the  determination  of  the  MMR.  The 
significant  wave  height  calculated  from  IRTOOL  was  1 .46  m  while  for  IRBLEM  it  was 
0.74  m.  When  we  edited  the  cut-off  height  in  IRBLEM  so  as  to  be  the  same  as  IRTOOL, 
IRBLEM  did  not  predict  the  existence  of  any  secondary  image  but  now  the  MIVRs  were 
changed  a  little  (Figure  67,  symbol 

The  infrared  (IR2)  camera  was  positioned  at  h,.=  7.3  m  above  the  water  level.  As 

we  can  see  from  Figure  70  only  one  MIVR  data  point  was  observed  for  the  target  at  14.5 
m  above  the  water  level  (ship’s  stack)  and  for  that  point  the  IRTOOL  predicted  a  MIVR 
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which  is  about  70  m  longer  while  IRBLEM  predicted  a  MIVR  800  m  shorter.  All  the 
profiles  for  IRTOOL  remain  the  same  as  before  while  for  IRBLEM  only  the  profiles  of 
the  refractivity  and  its  first  gradient  changed  (Figures  71 A  and  7 IB)  due  to  the  band 
dependent  coefficients  A,  B  that  are  used  for  the  calculation  of  refractivity.  As  we  can  see 
from  Figure  71 B  the  gradient  of  refractivity  in  IRBLEM  in  this  case  is  bigger  than  in 
IRTOOL.  That  means  the  ray  bending  is  greater  in  IRBLEM,  which  results  the  MIVRs 
and  MMRS  to  move  in  shorter  ranges.  Figure  72  shows  that  no  MMR  was  measured  for 
the  ship’s  stack  and  also  that  the  IRTOOL  did  not  predict  the  existence  of  any  mirages  for 
any  target  heights,  while  IRBLEM  predicted  the  existence  of  weak  secondary  images  for 
all  the  targets  above  about  4  m  from  the  water  level.  When  we  edited  the  cut-off  height  in 
IRBLEM  so  as  to  be  the  same  as  IRTOOL,  IRBLEM  did  not  predict  the  existence  of  any 
secondary  image,  but  the  MIVR  for  the  target  at  14.5  m  was  changed  slightly  (Figure  70, 
symbol 

Figure  48  shows  the  measured  ASTD  from  the  Hr.  Ms.  Tydeman  for  the  period 
from  03:10  to  05:25  GMT.  The  air  temperature  was  quite  steady  near  1 1.2 0  C  during  the 
track  while  the  sea  temperature  was  steady  near  1 1 .6 0  C  for  the  first  32  km  before 
jumping  up  to  12.4°  C.  As  no  lights  were  visible  beyond  32  km  the  slight  horizontal 
variation  in  ASTD  (0.4 0  C)  is  not  expected  to  have  a  significant  effect  on  the  agreement 
between  the  models  and  the  data. 
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Height  Above  WL  (m) 


Figure  22.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2010F2 
event  measured  by  the  visible  VI  camera.  The  Angular  Separation  is  the  elevation  angle 
of  each  image  relative  to  the  higher  target  source.  Each  image  in  the  upper  branch  is  a 
mirage. 
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Figure  23  A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 


Figure  23B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  23C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  23D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  23E.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Custom 
atmosphere)  and  IRBLEM  model. 
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Figure  23F.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  24.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  after  editing  the 
cut-off  height  of  IRBLEM,  for  the  P2010F2  event  measured  by  the  visible  VI  camera. 
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Figure  26.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2010F2 
event  measured  by  the  visible  VI  camera  for  ASTD  =  -3.9  °C. 
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Figure  27.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2010F2 
event  measured  by  the  visible  VI  camera  for  ASTD  =  -3.7  0  C. 
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Figure  29.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2010F3 
event  measured  by  the  visible  V 1  camera. 
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Figure  30A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 


Figure  30B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  30E.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  30F.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  31.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  after  editing  the 
cut-off  height  of  IRBLEM,  for  the  P2010F3  event  measured  by  the  visible  VI  camera. 
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Figure  32.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2810P2 
event  measured  by  the  visible  V2  camera. 


Figure  33A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 


Figure  33B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  33C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  33D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 


126 


P2810P2  (V2) 


- 

X 

X 

o  0 

1 

X 

o 

X 

- cr 

o 

\ 

£ 

X 

X 

o 

o 

rp 

4* 

£ 

®  12  60 

X 

X 

C 

C 

0 

x 

X 

X 

C 

c 

> 

3  1 

X 

< 

5  | 

O  IRTOOL  CL 

I3TOM  ATMOSPH 

ERE 

X 

V 

£  QQ  ___ 

X  IRTOOL  ML 

-  IRBLEM 

[w  ATMOSPHERE 

X 

X 

X' 

5 

D 

'“'i  £;-~ 

2.50  _ 

X 

X 

X 

l 

ss 

. .  XT 

J  o 

c .fj 

'J 

•T 

1 

1  1 

■  i 

277.00  278.00  279.00  280.00  281.00  282.00  283.00 


Refraetivity 


Figure  33E.  Refraetivity  profiles  as  predicted  from  IRTOOL  (using  the  Custom 
atmosphere)  and  IRBLEM  model. 
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Figure  33F.  Gradient  of  the  Refraetivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  34.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  after  editing  the 
cut-off  height  of  IRBLEM,  for  the  P2810P2  event  measured  by  the  visible  V2  camera. 
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Figure  36.  Comparison  of  the  LWKD  and  IRTOOL  model  MIVRs  with  the  ship-tracking 
data  measured  by  the  visible  (V2)  camera. 
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Figure  3  7  A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  37B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  37C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
and  Custom  atmosphere)  and  IRBLEM  model. 


Figure  37D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 


Figure  37E.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Custom 
atmosphere)  and  IRBLEM  model. 
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Figure  37F.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  38.  Comparison  of  the  LWKD  and  IRTOOL  model  MMRs  with  the  ship-tracking 
data  measured  by  the  visible  (V2)  camera. 
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Figure  39.  Comparison  of  the  LWKD  and  IRTOOL  model  MIVRs  with  the  ship-tracking 
data  measured  by  the  infrared  (IR1)  camera. 
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Figure  40A.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  40B.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  and  Custom  atmosphere)  and  IRBLEM  model. 
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Figure  41 .  Comparison  of  the  LWKD  and  IRTOOL  model  MMRs  with  the  ship-tracking 
data  measured  by  the  infrared  (IR1)  camera. 
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Figure  43  A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  43B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  43C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
atmosphere)  and  IRBLEM  model. 
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Figure  43D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  43E.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  43F.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 
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Figure  44.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  after  editing  the 
cut-off  height  of  IRBLEM  for  the  P2210J2  event  measured  by  the  visible  V2  camera. 
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Figure  46.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2510M1 
event  measured  by  the  visible  VI  camera. 


144 


Humidity  (%) 


Figure  47A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  47B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  47C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
atmosphere)  and  IRBLEM  model. 


Figure  47D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  47E.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  47F.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 
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Figure  48.  Variation  of  the  sea  and  air  temperature  along  the  ship’s  course  during  ship 
tracking  event  F. 
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Figure  50.  Comparison  of  IRTOOL  and  IRBLEM  model  calculations  for  the  P2810P3 
event  measured  by  the  visible  V2  camera. 
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Figure  51  A.  Humidity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  5 IB.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  5 IE.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  5  IF.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model.  . 
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Figure  53C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
atmosphere)  and  IRBLEM  model. 
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Figure  53D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  53E.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  53F.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 
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Figure  54.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship 
tracking  data  measured  by  the  infrared  (VI)  camera. 
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Figure  55.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship¬ 
tracking  data  measured  by  the  infrared  (IR1)  camera. 


Refractivity 

Figure  56A.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  56B.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 


Figure  57.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship¬ 
tracking  data  measured  by  the  infrared  (IR1)  camera. 
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Figure  58.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship 
tracking  data  measured  by  the  visible  (VI)  camera. 


Figure  59A.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
atmosphere)  and  IRBLEM  model. 
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Figure  59B.  Temperature  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  59C.  Pressure  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude  Winter 
atmosphere)  and  IRBLEM  model. 


Figure  59D.  Wind  speed  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  59E.  Refraetivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 
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Figure  59F.  Gradient  of  the  Refraetivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 
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Figure  60.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship 
tracking  data  measured  by  the  visible  (VI)  camera. 
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Figure  61 .  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship 
tracking  data  measured  by  the  visible  (V2)  camera. 
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Figure  62.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship 
tracking  data  measured  by  the  visible  (V2)  camera. 
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Figure  63.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship 
tracking  data  measured  by  the  visible  (V3)  camera. 
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Figure  65.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship 
tracking  data  measured  by  the  visible  (V4)  camera. 
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Figure  66.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship¬ 
tracking  data  measured  by  the  visible  (V4)  camera.  Neither  IRBLEM  nor  IRTOOL  shows 
a  mirage  for  this  case  even  though  we  had  one  observed  mirage  for  the  target  at  20.7  m. 
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Figure  67.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship¬ 
tracking  data  measured  by  the  infrared  (IR1)  camera. 
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Figure  68A.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


First  Derivative  of  Refractivity 


Figure  68B.  Gradient  of  the  Refractivity  profile  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 
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Figure  69.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship¬ 
tracking  data  measured  by  the  infrared  (IR1)  camera. 
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Figure  70.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MIVRs  with  the  ship 
tracking  data  measured  by  the  infrared  (IR2)  camera. 
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Figure  71  A.  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the  Mid-Latitude 
Winter  atmosphere)  and  IRBLEM  model. 


Figure  7  IB.  Gradient  of  the  Refractivity  profiles  as  predicted  from  IRTOOL  (using  the 
Mid-Latitude  Winter  atmosphere)  and  IRBLEM  model. 
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Figure  72.  Comparison  of  the  IRBLEM  and  IRTOOL  model  MMRs  with  the  ship¬ 
tracking  data  measured  by  the  infrared  (IR2)  camera.  IRTOOL  shows  no  mirage  for  this 
case  and  none  was  observed. 
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VIII.  CONCLUSIONS 


In  this  thesis  we  compared  the  refractivity  data  obtained  during  the  MAPTIP  trial 
at  Katwijk  Netherlands  in  1993  with  the  predictions  made  by  the  non-reff  active 
atmospheric  (geometric)  model  and  the  refraction  models  of  IRTOOL  and  IRBLEM.  Due 
to  the  large  amount  of  data  only  a  portion  of  these  data  with  the  concurrent 
meteorological  conditions  were  chosen  here,  so  as  to  cover  an  appreciable  range  of 
ASTDs  and  wind  velocities  sampled  at  various  times  of  the  day.  The  events  chosen  were 
often  observed  by  cameras  from  more  than  one  participating  nation. 

The  extended  research  showed  that  IRTOOL  could  predict  mirages  in  cases  in 
which  they  were  found  experimentally  if  the  code  cut-off  height  is  decreased 
significantly.  As  noted,  the  wave  height  is  an  important  factor  that  controls  the  Maximum 
Mirage  Height  (MXMH)  or  the  beginning  of  the  mirage  zone.  When  the  wave  height 
increases,  the  MXMH  decreases,  or  the  beginning  of  the  mirage  zone  moves  to  longer 
ranges  and  vice  versa.  In  the  ray-tracing  program  in  IRBLEM  the  cut-off  height  is 
controlled  by  the  Factor  Wave  Height  (FWH)  which  is  an  input  parameter.  The  value  of 
this  factor  can  be  very  useful  for  the  purposes  of  research  so  as  to  seek  the  ideal  cut-off 
height,  meaning  that  using  different  FWH  the  more  appropriate  cut-off  height  which  will 
correspond  to  the  real  one  can  be  estimated.  The  nominal  value  of  FWH  (0.6)  that  is  used 
in  this  thesis  results  in  a  cut-off  height  0.78HX ,  and  seems  to  be  very  accurate.  In  cases 
where  the  cut-off  height  used  by  IRTOOL  and  IRBLEM  was  almost  the  same  (events 
P2810P3  and  T2810),  the  results  given  by  the  two  models  were  almost  identical  and  in 
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very  good  agreement  with  the  observed  data.  In  all  other  cases  where  the  cut-off  height 
used  by  IRTOOL  was  much  greater  than  both  the  one  used  by  IRBLEM  and  the  measured 
value,  IRTOOL  seemed  to  give  poor  predictions  for  the  mirages.  The  equation  that 
IRTOOL  uses  to  calculate  the  cut-off  height  seems  to  be  inaccurate.  The  results  given  by 
IRTOOL  would  be  improved  significantly  if  this  equation  were  modified  or  even  better  if 
the  measurement  of  the  wave  height  was  an  input  in  the  IRTOOL  input  parameters  as  is 
done  in  IRBLEM.  A  deficiency  of  IRTOOL  is  that  in  some  cases  the  measurement  of  the 
wind  speed  is  taken  at  a  different  height  from  the  other  meteorological  measurements, 
while  the  program  assumes  that  all  the  measurements  were  taken  at  the  same  height.  This 
results  in  an  error  in  the  wind  speed  profile  and  consequently  in  the  cut-off  height,  the 
calculation  of  which  is  based  on  the  wind  speed  at  10  m  height. 

One  difference  between  the  two  models  is  that  the  pressure  is  not  an  input 
parameter  in  IRTOOL,  while  it  is  in  IRBLEM.  As  was  shown,  the  error  in  using  the 
midlatitude  winter  atmosphere  with  its  default  pressure  profile  is  no  greater  than  0.93% 
and  so  its  contribution  to  the  gradient  of  refractivity  is  neglible. 

Another  difference  between  the  two  models  is  that  IRTOOL  does  not  model 
dispersion  in  the  calculation  of  refractivity  as  it  only  uses  one  equation  for  all  the 
wavebands.  IRBLEM,  on  the  other  hand,  uses  four  different  equations  dependent  on  the 
wavelength.  From  these  four  equations  only  the  equation  that  is  used  for  the  far  infrared 
deviates  significantly  from  the  other  three  (the  equations  for  the  visible,  near  and  mid 
infrared)  and  from  the  equation  used  by  IRTOOL.  As  was  expected  for  the  visible,  near 
and  mid  infrared,  both  models  predicted  about  the  same  results  (for  the  same  cut-off 
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height),  which  were  in  very  good  agreement  with  the  observed  data.  Moreover,  in  the  far 
infrared  the  results  predicted  by  IRTOOL  were  much  better  than  those  predicted  by 
IRBLEM. 

The  general  impression  given  by  the  analysis  is  that  the  refractivity  profile,  which 
is  primarily  dependent  upon  the  temperature  profile,  produced  by  IRTOOL  causes  more 
ray  bending  than  that  produced  by  the  LWKD  model  of  IRBLEM.  Therefore  IRTOOL 
becomes  more  effective  in  the  prediction  of  mirages. 

From  the  results  of  the  comparison  with  the  MPN  observation  data,  both  models 
give  a  good  fit  to  the  experimental  data  under  the  assumption  that  both  use  the  same  cut¬ 
off  height  which  controls  the  MXMH.  Because  of  the  quality  of  the  fit  to  the  MPN  data, 
we  expect  the  fit  of  these  models  to  the  ship-tracking  events  to  be  reasonably  good.  In  the 
visible,  near  and  mid  infrared  the  MIVR  prediction  from  IRTOOL  was  always  within  2.7 
km  from  the  measured  data.  The  highest  deviation  in  the  results  occurred  when  the 
variation  in  the  ASTD  measured  by  the  ship  was  significant  (in  the  order  of  2 0  C).  In 
cases  where  this  uncertainty  was  small  (about  0.3 0  C)  the  MIVR  prediction  from  IRTOOL 
was  within  1.3  km  from  the  measured  data.  For  the  same  wavebands  the  MIVR  prediction 
from  IRBLEM  was  always  within  2.4  km  from  the  measured  data,  when  the  uncertainty 
in  ASTD  measured  from  the  ship  was  big  enough,  while  in  cases  where  this  horizontal 
inhomogeneity  was  small  the  MIVR  predictions  were  within  1  km  from  the  measured 
data.  In  the  far  infrared  IRTOOL  MIVR  predictions  were  much  better  than  IRBLEM.  In 
order  to  compare  the  MMR  results  we  make  the  assumption  that  both  models  use  the 
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same  cut-off  height,  namely  that  used  by  IRBLEM,  which  was  always  closer  to  the 
measured  HX  wave  height.  In  this  case  the  MMR  results  for  both  models  are  in  good 
agreement  among  themselves  and  with  the  experimental  data.  The  IRBLEM  results  are 
always  within  0.8  km  from  the  experimental  data  while  those  for  IRTOOL  are  always 
within  1 .2  km  and  more  often  than  not  well  within  0.5  km  of  the  measured  data.  These 
discrepancies  are  quite  reasonable  as  the  uncertainty  in  the  measured  range  was  about  200 
m.  The  ship  was  often  bouncing  up  and  down  in  the  waves  by  several  meters  so  that  for  a 
source  at  a  height  of  10  m  a  variation  in  height  of  1  m  could  lead  to  variation  in  the 
predicted  MIVR  of  about  500  m. 

Some  of  the  differences  between  the  predictions  obtained  from  the  different 
models  and  the  experimental  data  are  due  to  the  gradient  profiles  used  by  the  models  and 
particularly  the  temperature  profile  near  the  water  surface.  The  temperature  profile  could 
be  modified  so  as  to  better  describe  the  region  close  to  the  water  surface.  In  particular  a 
weaker  temperature  gradient  for  the  first  few  meters  would  allow  the  curves  to  move 
towards  the  measured  data.  Another  reason  that  we  have  already  mentioned  was  the 
horizontal  ASTD’s  inhomogeneity  measured  from  the  Hr.  Ms.  Tydeman.  Generally,  we 
can  say  that  variations  in  ASTD  affect  IRTOOL  more  than  IRBLEM.  Even  though  the 
results  from  IRBLEM  are  slightly  better  than  IRTOOL,  we  must  note  that  in  most  of  the 
cases  the  humidity  calculated  from  IRBLEM  at  the  sea  surface  was  greater  than  100% 
which  is  a  non-physical  result.  Also,  the  water  temperature  used  by  IRBLEM  does  not 
correspond  to  the  measured  value  even  though  the  measured  value  is  an  input  in  the 
IRBLEM  parameters. 
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The  comparison  of  the  geometric  models  of  the  ship-tracking  events  between 
IRBLEM  and  IRTOOL  showed  that  in  all  cases  IRBLEM  calculates  a  geometrical  range 
greater  than  IRTOOL.  We  showed  that  this  happens  because  IRBLEM  does  not  take  into 
account  the  wave  height  for  the  calculation  of  the  geometrical  optical  sight  range. 

Finally,  it  is  reasonable  to  conjecture  whether  the  differences  in  the  results 
produced  from  the  ray  tracing  programs  of  IRTOOL  and  IRBLEM  are  dependent  only 
upon  the  refractivity  profiles.  Another  factor  that  might  produce  significant  discrepancies 
between  the  predictions  of  the  two  models  is  the  difference  in  the  ray  tracing  algorithms. 
To  remove  such  doubts,  several  refractivity  profiles  produced  by  the  two  models  for  a 
number  of  cases  should  be  used  as  an  input  to  a  common  ray-tracing  program  and  their 
corresponding  outputs  should  be  compared.  This  could  probably  be  the  objective  of  a 
future  research.  However,  in  this  thesis  we  showed  that  in  cases  such  as  the  events 
P2810P3  and  T2810K,  where  the  profiles  of  the  gradient  of  refractivity  were  very  similar 
and  also  the  cut-off  height  that  was  used  from  each  model  was  almost  the  same,  the 
results  predicted  from  the  two  models  were  in  extremely  good  agreement.  So,  it  is  safe  to 
say  that  the  differences  in  the  results  produced  by  the  two  models  appear  to  be  solely  due 
to  the  differences  between  the  profiles  generated  from  IRTOOL  and  IRBLEM. 
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